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Abstract 

We review the theoretical motivations and experimental status of searches for stable mas- 
sive particles (SMPs) which could be sufficiently long-lived as to be directly detected at 
collider experiments. The discovery of such particles would address a number of impor- 
tant questions in modern physics including the origin and composition of dark matter in 
the universe and the unification of the fundamental forces. This review describes the tech- 
niques used in SMP-searches at collider experiments and the limits so far obtained on the 
production of SMPs which possess various colour, electric and magnetic charge quantum 
numbers. We also describe theoretical scenarios which predict SMPs and the phenomenol- 
ogy needed to model their production at colliders and interactions with matter. In addition, 
the interplay between collider searches and open questions in cosmology is addressed. 
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1 Introduction 



An open question in modern physics is the possible existence of heavy, exotic par- 
ticles, which can be directly detected at collider experiments through their interac- 
tions with matter. This paper describes collider searches for these so-called Stable 
Massive Particles (SMPs), together with a review of theoretical models in which 
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such states appear. We also discuss the astrophysical consequences that an SMP 
discovery would imply. In this work, we define an SMP as a particle which does 
not decay during its passage through a detector, and which would undergo electro- 
magnetic and/or strong interactions with matter. Normally we would expect such 
particles to be heavier than a proton. 

Historically, the strange, long-lived kaons heralded a revolution in particle physics, 
in terms of new fundamental matter, and also in the shape of a new (approximately) 
conserved quantum number. Today, states classifiable as SMPs recur in many the- 
oretical extensions of the Standard Model (SM). One crucial aspect in this context 
is naturally the question of dark matter in the universe, but also charge quantisa- 
tion (magnetic monopoles? millicharges? grand unification?), the flavour question 
(more generations?), parity violation (vector-like fermions? mirror fermions?), the 
hierarchy problem (supersymmetric SMPs?), and more (SMPs from extra dimen- 
sions? your favourite particle?) are relevant. We devote Section 2 to a review of 
theoretical scenarios with SMP states. 

Due to the symbiosis between accelerator particle physics and astrophysics for 
long-lived particles, we focus on the cosmological and astrophysical implications 
of SMPs in Section 3. This discussion includes the implications of an SMP discov- 
ery at an accelerator on topics such as dark matter and nucleosynthesis. 

In most non-generic searches for SMPs, theoretical models are needed to calculate 
their production cross sections and final-state topologies. In the case of massive 
coloured objects, their fragmentation into jets must be described. Section 4 contains 
an overview of the techniques used to model the production of SMPs. 

The detection of SMPs is only possible once the interactions in matter are under- 
stood. In Section 5 we summarise the phenomenology used to model the interac- 
tions of SMPs with detector material, including a description of their electromag- 
netic and hadronic interactions. 

In Section 6, a range of possible search techniques which can be used to iden- 
tify SMPs are described. These include techniques based on ionisation energy 
losses, Cherenkov techniques and time-of-flight methods. Additionally, specific 
techniques used in magnetic monopole searches are summarised. 

A wide variety of SMP searches have been performed in lepton-hadron, electron- 
positron and hadron-hadron reactions [1], yielding important constraints on the pa- 
rameter spaces of a number of theories. The searches vary from the very general 
types, which make minimal assumptions regarding the quantum numbers of the 
SMP, to those performed within a specific theoretical scenario. Section 7 presents a 
summary of the most important model-independent searches, as well as a selection 
of searches made within commonly studied exotic physics scenarios. 

Approaching the exciting time when the Large Hadron Collider (LHC) will pro- 
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duce its first collisions, it is appropiate to discuss prospects for the discovery of 
SMPs with this machine. Searches for SMPs with masses up to several TeV will be 
possible, representing an order of magnitude increase in mass sensitivity compared 
with earlier colliders. Section 8 discusses the the discovery potential of a range of 
types of SMPs at the LHC. Finally, Section 9 contains a summary of this report. 

Our report is complementary to, and differs from, previous reviews of SMPs [2-5] 
in a number of ways. The principal difference is that we focus on collider searches, 
although non-accelerator results such as those obtained using mass spectroscopy 
and cosmic ray detection are discussed, when appropriate. In view of the circa 50 
accelerator searches, which have taken place over the past two decades, an up- 
to-date and detailed description of the techniques and results of these searches, 
and of their theoretical motivation, is overdue. By specialising in collider studies 
we are able to cover searches for many different proposed species of SMPs and 
the experimental challenges specific to each of them. For example, we provide a 
detailed treatment of hadronic SMPs, which is lacking in earlier reviews. Further- 
more, although magnetic monopoles are usually treated separately to electrically 
charged SMPs in this type of article, they are included here since both types of 
particles are expected to share many common experimental signatures. Magnetic 
monopole searches can provide sensitivity to SMPs with values of electric charge 
beyond those covered by dedicated searches for electrically charged SMPs, and 
vice versa. Thus, by considering a broader range of collider searches, it is possible 
to more fully describe the types of SMPs, and their possible masses, which have 
been excluded, and those which could have been produced at colliders but which 
potentially remain unobserved. 



2 Theoretical scenarios for SMPs 

In this section, we first give a brief introduction to theoretical possibilities for SMP 
states, followed by in-depth discussions on some of the more commonly considered 
scenarios. In particular we discuss supersymmetric models in Section 2. 1 and mod- 
els with universal extra dimensions in Section 2.2. Sections 2.3 and 2.4 concern 
alternative possibilities that do not fit smoothly into either of the above categories, 
such as -induced millicharges and magnetic monopoles. 

To gain an impression of the wide range of possible SMP states, it is instructive first 
to take a look at the Standard Model itself. Protons, neutrons, electrons, and muons 
are all examples of stable particles which undergo interactions in a detector. So are 
TT"*", K'^, and K^. These well-studied states also illustrate very well the spectrum 
of possibilities beyond the Standard Model. Consider the following well-known 
properties, and apply them implicitly to hypothetical new states. The electron is a 
fundamental particle. It does not decay, e.g. to neutrinos due to the conservation of a 
gauged quantum number, electric charge in this case. It is the lightest state carrying 
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that quantum number. The proton is a complicated bound state composed of a set 
of more fundamental particles which are held together by a high- strength short- 
range gauge force. The kinematically allowed decay of the proton to, e.g. e+7r° 
does not occur due to conservation of certain global (i.e. non-gauged) quantum 
numbers, in this case baryon and lepton number. The conservation of these numbers 
are accidental symmetries of the model. The long lifetime of the neutron owes to 
a combination of the weak force being involved and the very small decay phase 
space. Though unstable in free space, there are still many of them around, due to 
the existence of stable "composites of composites", the nuclei of baryonic matter. 
The muon has a comparatively long lifetime due to the hierarchy between the muon 
mass and the weak scale (F^ oc m^/m^). At low energies, the decay appears to 
proceed via a non-renormalisable dimension-6 operator which, in the fundamental 
theory, arises due to a virtual massive gauge boson, the W. Finally, also the atom 
furnishes an example. Its stability could not be understood in terms of Maxwell's 
theory. Only after a drastic revision, quantum mechanics, could it be accounted for. 
The corresponding case today would be the discovery of a state whose stability 
could not be accounted for within the framework of quantum field theory. 

Thus, already the SM contains quite a varied history of stable interacting states. 
Turning now to physics beyond the SM, there are several generic possibilities for 
SMP's, essentially all recurrences of the states mentioned above, with the addition 
of topological defects, like cosmic strings or magnetic monopoles. 

The most obvious possibility for an SMP is that one or more new states exist which 
carry a new conserved, or almost conserved, global quantum numbefl SUSY with 
i?-parity, extra dimensions with KK-parity, and several other models fall into this 
category. The lightest of the new states will be stable, due to the conservation of 
this new "parity", and depending on quantum numbers, mass spectra, and interac- 
tion strengths, one or more higher-lying states may also be stable or meta-stable. 
In general, electrically charged stable states are excluded by cosmology, and also 
coloured particles are strongly constrained, as will be discussed in Section 3. For 
this reason, and to obtain a solution to the dark-matter problem, models are usually 
constructed to provide Mn-charged stable dark-matter candidates, most often in the 
guise of weakly interacting massive particles (WIMPs). From a motivational point 
of view, SMP models thus come in two categories: 

(1) Models which solve the dark-matter problem with a WIMP-type dark-matter 
particle, but which also have one or more higher-lying meta-stable SMP states. 

(2) Models which have long-lived SMP states, but which either do not address 
dark matter or address it with a non-WIMP dark-matter particle. 



As an aside, at the most fundamental level any such global symmetry probably has to 
be the remnant of a broken gauge or space-time symmetry (e.g. KK parity is the remnant 
of broken higher-dimensional Poincare invariance), to avoid stability problems [9-11], but 
we shall here treat them simply as discrete global symmetries, whatever their origin. 
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As a quick reference, Tabs. 1 and 2 give a condensed overview of the SMP states 
discussed in the text, along with a description of which scenario(s) give rise to each 
of them. 



2.1 SMP states in Super symmetry 

Among the most interesting and comprehensively studied possibilities for observ- 
able new physics is supersymmetry (SUSY) — for excellent reviews see, e.g. [12, 
13]. Interesting in its own right by being the largest possible space-time symmetry, 
it was chiefly with a number of additional theoretical successes during the eighties 
that supersymmetry gained widespread acceptance, among these successful gauge 
coupling unification, a natural candidate for dark matter, and an elegant solution to 
the so-called hierarchy problem associated with the smallness of Mz/M^-^^^^y.. 

Stated briefly, supersymmetry promotes all the fields of the SM to superfields. Each 
{N = 1) superfield contains one boson and one fermion as dynamical degrees of 
freedom. Hence not only should all the SM particles (including at least one more 
Higgs doublet) have 'sparticle' superpartners, with spins differing by 1/2, but also 
the interactions of these superpartners should to a large extent be fixed by super- 
symmetry. However, due to the absence of observed sparticles so far, there must 
exist a mass gap between the SM particles and their superpartners, implying that 
supersymmetry must be broken at the weak scale. This mass gap cannot be much 
larger than ~ 1 TeV without invalidating the solution to the hierarchy problem. 
Independently, both SUSY dark matter and SUSY grand unification also point to 
roughly the TeV scale, hence there is good reason to believe that SUSY, if it exists, 
should manifest itself at the next generation of colliders, and a long list of sophisti- 
cated theory tools for SUSY phenomenology are available, see, e.g. [14-18]. 

In the most general SUSY models, both lepton and baryon number are violated. 
If the SUSY-breaking mass scale [19] is 0{\ TeV), then bounds on proton decay 
[1,20] require the product of the B and L violating couplings to be less than 10^^^ 
[21]. Hence, either one or both of the accidental symmetries that protect B and 
L in the Standard Model are usually promoted to explicit ones in SUSY models. 
The most common way to accomplish this is through postulating the conservation 
of the quantum number R = (_i)3B+l+25 [22], where B(L) is baryon (lepton) 
number and S is spin. Since SM particles carry R = +1 and their SUSY partners 
R— —1, the lightest supersymmetric particle (LSP) is stable, providing a possible 
solution to the dark-matter problem. If i?-parity is violated, either lepton or baryon 
number must still be conserved, to avoid rapid proton decay. In this case, the LSP 
is not stable (though it may still be long-lived) and the dark-matter problem is not 
addressed. For a charged or coloured LSP, i?-parity violation can be invoked to 
obtain a lifetime in the allowed range [7, 21]. 
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SMP LSP Scenario Conditions 





X? 


MSSM 


fi mass (determined by m?^ ^, /x, tan and A^) close to Xi 
mass. 




G 


GMSB 
^MSB 
SUGRA 


Large N, small M, and/or large tan j3. 

No detailed phenomenology studies, see [23]. 

Supergravity with a gravitino LSP, see [24]. 




fi 


MSSM 
AMSB 

5MSB 


Small rufj^ ^ and/or large tan /? and/or very large At-. 
Small mo, large tan/3. 
Generic in minimal models. 


In 


G 


GMSB 


fi NLSP (see above), ei and /ii co-NLSP and also SMP for 
small tan /3 and jj,. 






5MSB 


(=1 and jli co-LSP and also SMP when stau mixing small. 


xt 


x2 


MSSM 
AMSB 


m ~+ — m^o ^ m^+. Very large Mi ^2 ~ 2 TeV ^ (Hig- 
gsino region) or non-universal gaugino masses Mi k, AM2, 
witn tne latter condition reiaxea to ivii <, IV12 lor IVI2 
Natural in 0-11 models, where simultaneously also the g can 
be long-hved near 6gs = —3. 

Ml > M2 natural, mo not too small. See MSSM above. 
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X? 


MSSM 


Very large m| » M3, e.g. spht SUSY. 




G 


GMSB 


SUSY GUT extensions [25-27]. 
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MSSM 
GMSB 


Very small M3 < Mi,2, O-II models near (5gs = -3. 
SUSY GUT extensions [25-29]. 


ii 


X? 


MSSM 


Non-universal squark and gaugino masses. Small m? and 
M3, small tan /3, large At. 


hi 






Small m| and M3, large tan P and/or large A), 3> At. 



Table 1 



Brief overview of possible SUSY SMP states considered in the literature. Classified by 
SMP, LSP, scenario, and typical conditions for this case to materialise in the given scenario. 
See text for details. 

2.1.1 The MSSM and CMSSM 

At present, little is known for certain about the nature of the SUSY-breaking mech- 
anism. Some special cases will be discussed below, but at the most general level we 
define the minimal supersymmetric Standard Model (MSSM) to contain all pos- 
sible SUSY-breaking interactions, which are consistent with gauge and Poincare 
invariance, and which do not cause the hierarchy problem to reappear (so-called 
soft SUSY-breaking). This implies an additional mass term for each of the gaugi- 
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nos Mi,2,3, extra mass terms squared mfj for each pair of scalars i and j, and also 
trilinear couplings Aij^ between combinations of three scalars i, j, and A;. Due to 
significant constraints on CP and flavour violation, these parameters are normally 
taken to be real and flavour diagonal. To further limit the number of parameters, 
the additional assumption of GUT scale unification of the scalar masses mo, gaug- 
ino masses mi/2, and trilinear parameters Aq is often made, a choice sometimes 
referred to as the constrained MSSM or CMSSM. In addition, there are also the 
SUSY-conserving Yukawa couplings Yij^, gauge couplings gi, vacuum expecta- 
tion values (vev's) for each of the Higgs doublets vi^2, the n parameter (a SUSY- 
conserving mixing term between the two Higgs doublets), and, if i?-parity is vi- 
olated, Yukawa-like interactions between lepton and quark superfields as well as 
L-violating mixing terms between one of the Higgs superfields and the lepton su- 
perfields. The sum of the Higgs vev's is fixed by the Z mass, and their relative size 
is usually cast in terms of the parameter tan/3 = V2/V1, with tan/3<50 required 
for perturbativity up to the GUT scale. See [12, 13] for further details. 

As Tab. 1 suggests, the general MSSM allows for essentially any sparticle to be 
an SMP. In models with neutralino dark matter, the NLSP can be long-lived if its 
decay phase space is small or zero, mimicking the case of the neutron in the SM. An 
interesting scenario here is that of a light stop NLSP, as motivated by electroweak 
baryogenesis [30, 31]. In this scenario non-universal squark mass terms are used to 
arrange a small mass difference between the ii and the LSP Xi, while the lightest 
chargino is kept too heavy for the decay ii — > bxt to occur. In this case, only the 
radiative process ii ciXi is open, and the ii can be quite long-lived. See [30-33] 
for detailed studies. 

Though not equally well motivated at present, a light NLSP fi (small slepton 
masses and/or large tan l3 and/or large Ar) or less likely 61 (small squark masses 
and large tan /3 and/or very large Ah 3> At) could also be possible. For a fi NLSP 
with small mixing (small tan /3 and A^.), the ei and fii could simultaneously also 
be long-lived. A small mass difference m -+ — m^o < 171^^+ can occur for very large 
Mi 2~2 TeV ^ (Higgsino region) or non-universal gaugino masses [34-36] 
Ml > 4M2, with the latter condition relaxed to Mi > M2 for M2 < To obtain a 
long-lived gluino NLSP, the requirement is small M3, possibly in combination with 
very large m?, since the gluino decay proceeds via intermediate squarks (see, e.g. 
the discussion of split supersymmetry below). 

If the LSP is not required to be a neutralino, essentially all of the cases just men- 
tioned have obvious extensions where the NLSP becomes the LSP. In this case, 
dark matter is much more problematic, and i?-parity violation may be necessary to 
render the LSP unstable. 

When gravity is included, as in supergravity models, also the gravitino can in prin- 
ciple be the LSP, depending on its mass. If the gravitino is light, the NLSP will 
then often be long-lived due to the smallness of the gravitational coupling. For su- 
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pergravity models with a long-lived f NLSP see [24] ; other models with a gravitino 
LSP will be discussed below. Finally, we note that a similar phenomenology can 
also be obtained in extensions with an axino LSP [6]. 



2.1.2 Gauge Mediated Supersymmetry Breaking 

In models with gauge mediation [37, 38] (GMSB, see [39^1] for mass spectra 
and [42] for a review), the gravitino is very light (m^ ^ 1 keV) and hence the LSP 
for any relevant choice of parameters. Minimal models are cast in terms of six pa- 
rameters, typically A^, A, M, tan/?, sgn(yu), and Cgrav Briefly stated, the meaning 
of these parameters is that N chiral SU (5) multiplets are added to the theory at the 
scale M ^ 10^^ GeV [42]. If not only SU (5) multiplets are added, the counting gets 
more complicated, but there is still an effective N. These "messengers" couple di- 
rectly both to the MSSM fields (via the ordinary SM gauge interactions) and also to 
an unspecified source of SUSY-breaking. A ~ 10 — 100 TeV is the effective SUSY- 
breaking scale, related to the fundamental SUSY-breaking scale ^/F by a relation 
A = F/M. The next-to-lightest sparticle (NLSP) decays only via the gravitational 
coupling and can be very long-lived. For a slepton NLSP [42]: 



Since the SUSY-breaking terms are induced by gauge interactions, they are flavour 
universal and their sizes are proportional to the amount of gauge charge carried 
by each field. In minimal models, the next-to-lightest supersymmetric particle is 
therefore always the "least charged" of either the gauginos or the scalars. The lat- 
ter occurs in particular for large values of the messenger index N . In a tiny and 
near-excluded parameter region at small values of the model parameters A and M, 
the NLSP is then a sneutrino, otherwise it is the fi (for a benchmark and model 
line, see Snowmass point 7 [43]). Note, however, that N cannot be chosen arbi- 
trarily large, perturbativity of the theory up to the GUT scale requiring, e.g. ^ 5 
for smaU M < 10^ GeV and A^ < 10 for M - 10^° GeV. To iUustrate the parameter 
space. Fig. 1 shows the smallest messenger index A^ required to have a fi NLSP 
for > and a) relatively light messengers M = 2A and b) heavy messengers 
M — 10^° GeV, as a function of tan P and A (since Cgrav only affects the decay 
we leave it unspecified). The numbers were obtained with ISAJET v.7.71 [44] using 
m,t = 175 GeV. The light grey areas at small A are theoretically excluded due to un- 
stable vacua and/or non-perturbative couplings at the GUT scale; no experimental 
or indirect constraints were included here. 




mm 



(1) 



with the gravitino mass controlled by Cgrav < 1 and F: 




(2) 
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A [TeV] 

b) 

Fig. 1. GMSB: the smallest index number N required to obtain a fi NLSP as a function of 
A and tan /? for aj light messengers (M = 2 A) and heavy messengers (M = 10^'^ GeV). 
The colour coding is the same for both plots and corresponds to the legend shown with b). 

With regard to dark matter, it is interesting to note that even in the presence of B- 
violation, the gravitino would still be stable and a dark-matter candidate, since no 
kinematically allowed decays would be available. For other possibilities for GMSB 
dark matter, see [42]. 

Finally, a long-lived fi is not the only SMP possibility in GMSB. If the mix- 
ing and consequently the mass splitting in the stau sector is not too large (small 
tan [45]), then the ei and jli may be nearly mass-degenerate (co-NLSP) with 
the fi and hence can simultaneously be SMPs. As with all supersymmetric scenar- 
ios, there is also a very large space of possible non-minimal models, cf. [42]. Of 
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m (GeV) 



Fig. 2. The gluino lifetime in split SUSY as a function of the scalar mass parameter m 
for tan/3 = 2, /i > 0, and various choices of the gluino mass, as calculated by [48]. The 
dashed horizontal line indicates the age of the Universe, tu = 14 Gyr. 

particular interest here are SUSY GUT extensions of GMSB in which the coloured 
messengers are naturally much heavier than their weak counterparts, resulting in a 
gluino NLSP [25-27] or even LSP [25-29], depending on the gravitino mass. 



2.1.3 Split Supersymmetry 

IMeta-stable coloured sparticles also arise in the so-called split SUSY scenario 
[46,47], in which all the scalars (except the ordinary Higgs) have very large masses, 
while the gaugino and higgsino masses remain at or around the weak scale. Though 
the hierarchy problem is not addressed (except anthropically), this naturally sup- 
presses both proton decay and CP and flavour violation. Since gluinos can only 
decay via squarks (independently of whether _R-parity is conserved or not), the 
gluino lifetime can be very large in this scenario, somewhat similar to the case of 
the muon in the SM. The competing channels are tree-level 3-body decays to two 
quarks plus a chargino or neutralino and radiative 2-body decays into a gluon plus 
a neutralino, see [48] for explicit calculations or [49] for a simplified treatment. For 
illustration, in Fig. 2 we include a plot from [48] showing the gluino lifetime as a 
function of the scalar mass parameter m for tan P = 2, fi > 0, and various choices 
of the gluino mass. 

Another interesting possibility in split SUSY, which is not ordinarily viable, is that 
SUSY-breaking could be communicated directly at tree-level [50, 51]. In this case, 
the scalar masses m could be very close to the fundamental SUSY-breaking scale 
\/F. Since the gluino decays via virtual squarks discussed above are suppressed by 



10 



171 ^ and the coupling to the gravitino goes like l/F, this would open the possibility 
for a large branching fraction for gluinos to gravitinos, g — > gG [48]. 



2.1.4 Other SUSY-breaking scenarios: AMSB, O-II, andgMSB 

Anomaly mediated SUSY-breaking (AMSB, see [52-54]) is a variant of supergrav- 
ity, where the explicit SUSY-breaking terms are 'switched off (or at least heavily 
suppressed), leaving a scaling (conformal) anomaly in the supergravity Lagrangian 
as the sole source of supersymmetry breaking. This "anomaly-mediated" contribu- 
tion is always present, but is usually much smaller than other supergravity terms. 
Akin to GMSB, its virtue is that it is flavour universal by construction. In its pure 
form, however, AMSB gives rise to tachyonic sleptons (negative m|), hence, to ob- 
tain a viable phenomenology, additional positive lepton mass squared terms have to 
be included. The parameters of a minimal model are thus the gravitino mass 7713/2, 
the usual tan/3 and sgn(//), as well as additional soft SUSY-breaking masses for 
all sfermions 7770- The ratio of gaugino masses at the weak scale is approximately 
I Ml I : IM2I : IM3I ~ 2.8 : 1 : 8 [54], making the lightest chargino and neutralino 
nearly mass-degenerate and wino-like. For large 7770, the neutralino is the LSP, and 
hence the chargino can be long-lived. For small 777o the LSP can be the fi, except at 
small tan (3 and 7773/2 where it is the Ufi [54]. For an LHC phenomenology study, 
see [55]. 

Another interesting explicit realisation of non-universal gaugino masses is fur- 
nished by the so-called O-II orbifold model, a string-inspired scenario in which 
supersymmetry breaking is dominated by the overall 'size' modulus field arising 
from the orbifold compactification [56,57]. Gaugino masses arise at one loop, and 
a large degree of non-universality is generated naturally. The free parameters of 
a minimal model are 7770, tan/3, sgn(/x), and 6gs- The latter is called the Green- 
Schwarz mixing parameter and preferentially lies in the range —5^5gs < 0, with 
the negative integers -4 and -5 preferred for the model studied in [57]. The mass 
spectrum depends sensitively on Sqs but a typical feature, unless |5gs| is very large, 
is I Ml I ^ IM2I, resulting in a near-degeneracy between the lightest neutralino and 
chargino, both of which will be wino-like as in AMSB. Consequently, the chargino 
can be very long-lived. For Sqs — —3, the gluino mass is zero at the high scale, 
M3 = 0, and hence close to this value the gluino is typically the LSP [26, 57] 
or, slightly farther away from the minimum, nearly mass-degenerate with both the 
Xi and the xl- Note, however, that the other gauginos also have minima around 
Sqs = —3, so the entire gaugino spectrum can become very light in this region. 

In gaugino-mediated SUSY-breaking [58,59] (^MSB), the MSSM is embedded in a 

5-dimensional compactified braneworld setup. The only non-vanishing soft SUSY- 
breaking terms at the compactification scale, Mc, are the gaugino masses. In the 
original models, Mc = Mqvt, which gives a fi LSP. The jli and ei are typically 
slightly heavier, due to the smaller mixing in the first two generations, but they 
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can also be SMP candidates. Due to the obvious constraints from cosmology, most 
subsequent attempts have focussed on finding models beyond the minimal where a 
WIMP-like LSP could be recovered. In particular, if > Mgut> then the extra 
GUT running from M^. (where the boundary conditions are true) to Mqut (below 
which the SM running takes over) can generate non-zero scalar masses at the GUT 
scale [60-62], which may push the fi mass above the x\ one. Another way of 
modifying the spectrum is by introducing non-universal gaugino masses, e.g. as 
in higher-dimensional GUT models where the extra dimension is larger than the 
inverse GUT scale, < Mqut [63], but except for small parameter regions the 
fi LSP still dominates. A Tevatron study of gaugino mediation can be found in [64]. 
The interesting case of a GMSB-like phenomenology could also be possible, with a 
gravitino LSP and the fi the NLSP [23], but detailed phenomenology studies have 
so far not been carried out. 



2.2 SMP States in Universal Extra Dimensions 

In models of Universal Extra Dimensions (UED) all the fields of the SM, including 
both matter and forces, are allowed to propagate in some number of extra dimen- 
sions, usually taken to be one or two. These models [65-67] provide an interesting 
scenario for TeV scale physics and are consistent with low-energy constraints [68]. 
The compactification is constructed such that momentum conservation in the ex- 
tra dimensions is preserved at tree level, leading to a discrete Kaluza-Klein (KK) 
quantum number in the effective 4D theory. This quantum number is then broken 
to a KK parity at the loop level [67]. At least the lightest of the KK excitations of 
SM particles will thus be stable and can be a dark-matter candidate [69-71]. 

In the simplest analysis, all the KK modes of the light SM particles (the photon, 
gluon, and first generation fermions) may be sufficiently long-lived to be SMP can- 
didates as well (they are stable at tree level and are almost mass-degenerate level 
by level), but note that the presence of non-zero boundary terms in the extra di- 
mension(s) can change this picture drastically [72]. If small, these boundary terms 
can be treated perturbatively [73], but they are not calculable from fundamental 
principles and should be considered free parameters of the theory. 

The Minimal UED (MUED) model with one extra dimension can be specified in 
terms of three parameters, R, A, and niH, where R is the size of the extra dimension 
(an /Z2 orbifold stretched from to vri?), A > is a scale at which the 
boundary terms mentioned above are assumed to vanish, and mn is the mass of 
the (ground state) Higgs boson. The viable range of R found in [65] is 300 GeV < 
R~^ < 10 TeV. For two extra dimensions [74] the lower bound becomes more 
dependent on the compactification and ranges between 300 and 700 GeV, while for 
more extra dimensions a bound cannot be reliably estimated. See [75] for a Monte 
Carlo implementation of this model. 
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Finally, we note that a GMSB-like phenomenology with a stable KK graviton and 
a meta-stable next-to-lightest KK particle, e.g. ri, has also been considered, see 
[71,76]. 

2.3 SMP States in other scenarios beyond the Standard Model (BSM) 

We now turn to a brief overview of more exotic possibilities for SMP states in BSM 
physics. In particular, it is interesting to investigate the extent to which other models 
which have been proposed to address the dark-matter problem can also give rise to 
SMP states. We do not claim to be complete — model space is in principle infinite 
— nor do we give a detailed discussion of each scenario, but we hope to illustrate 
the spectrum of ideas, and where possible point the reader to relevant literature 
where further information can be found. Also keep in mind that for essentially all 
these scenarios, with the exception of leptoquarks and to a lesser extent Little Higgs 
models, no convenient package of collider phenomenology tools yet exists. 

2.3.1 Models with parity-like symmetries 

A class of models which have recently attracted attention is Little Higgs with T- 
parity (for recent reviews, see, e.g. [77, 78]). As with several of the other proposals 

discussed here, T-parity (not to be confused with time reversal) serves the dual pur- 
pose of simultaneously suppressing contributions to electroweak precision observ- 
ables [79] and providing a WIMP-like dark- matter candidate [80], the "Lightest T- 
odd particle" (LTP). In the context of SMPs, however, minimal models do not have 
much to offer. The only comer of parameter space where the LTP can be charged 
and/or coloured [80], at large values of the symmetry-breaking scale / and small 
T-odd fermion masses fh, roughly / ^ 1 TeV and m ^ 300 GeV, is excluded from 
Tevatron squark searches. It should be possible to construct more viable models 
with SMPs using non-universal rficoioured > "^ieptons> but so far no explicit models 
in this direction have been constructed. 

More exotic possibilities for supersymmetric SMP candidates also exist, in partic- 
ular in models where the MSSM gauge groups and particle content are extended to 
include more superfields. A recent example is a variant of the so-called Fat Higgs 
model, based on the MSSM with an extra confining SU{3) symmetry [81], the "Fat 
Higgs with a Fat Top" [82]. In this model, quasi-stable exotic chiral superfields 
(i.e. a complex scalar and a fermion) appear, which are charged under a global Z2 
symmetry, which makes them approximately stable. The strongly interacting ones 
are probably outside the range of colliders, but there is a weak-scale electrically 
charged multiplet whose members have SMP properties [82]. 

In warped extra dimensions with "GUT parity" [83-86], the combination of extra 
dimensions and an effective TeV scale supersymmetric grand unification results in 
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KK towers not only of the SM gauge and Higgs fields but also of their SUSY-GUT 
partners, including XY bosons and coloured Higgs multiplets. In models where the 
SM fermions also propagate in the bulk, add KK-towers for these and their SUSY- 
GUT partners as well. A parity can be chosen such that the MSSM particles are even 
and their GUT partners odd, hence the lightest "GUT-odd" particle (LGP) is stable 
or long-lived if this quantum number is conserved or approximately conserved, 
respectively. In the earliest scenario of [83], the LGP is typically a light isospin-up 
(-down) colour triplet XY gaugino, with electric charge -1/3 (-4/3), but in the more 
recent models [84-86] a wide range of possibilities are open. 

In a recent 5D model [87] of dynamical SUSY-breaking (DSB), TeV scale exotic 
scalars with the quantum numbers of GUT XY bosons appear, so-called xyons. If 
a condition similar to ii-parity holds in the DSB sector, these states are long-lived. 
Their precise quantum numbers depend on the details of the DSB — in general 
they are both coloured and charged. In the simplest SU{5) case they lie in a colour 
triplet isospin doublet with electric charges Q = —1/3,-4/3, but also 5*0 (10) 
assignments are possible, for example an additional doublet with Q = 1/3,-2/3, 
could easily be possible. See [87] for mass spectra and phenomenology. 

Finally, there also exist a few more general ideas for possible (quasi-)stable BSM 

particles, including long-lived leptoquarks (see, e.g. [88]) and additional (gener- 
ations of) fermions [89, 90]. The latter can either be straightforward additions to 
the SM generations, e.g. a 4th generation with 4th flavour approximately conserved 
[91], or they can have a non-SM like structure. Mirror fermions (see, e.g. [8,92,93]) 
are extra fermions whose right-chiral members lie in SU{2) doublets while the 
left-handed ones are singlets, i.e. opposite to the SM. A 'vector-like' generation is 
comprised of an extra SM generation together with its mirror, as e.g. in = 2 
supersymmetric models [92,94]. 



2.4 Models with gauged symmetries 

When moving from global symmetries to local ones, several new possibilities open 
up. Even though explicit models are scarce, we shall still try to provide a reasonable 
overview and discussion. 

Consider first the case of a new U{l)x gauge group. If X is conserved then a cor- 
responding new massless exotic photon exists, 7' (or "paraphoton" [95]), which 
mediates a new long-range gauge force between particles charged under X. If any 
SM field is charged under X, then X must he B — L [96], and the coupling is con- 
strained to be extremely tiny gs-L ~ 10^^^ [97]. Furthermore, only a state carrying a 
non-SM-like combination of quantum numbers could be stable, since it would oth- 
erwise decay to SM particles via prompt photon (or W/Z) emission. The stability 
would hence not rely directly on the smallness of Qb-l but rather on the conser- 
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Examples of possible SMP states in a variety of models beyond the MSSM (for MSSM 
SMPs, see Tab. 1). Classified by electric charge Q, colour representation CQCD, spin S, 
and scenario. 



vation of other quantum numbers. We shall therefore not dwell on this possibility 
further. 

The simplest is thus to postulate a genuinely new U{l)x, with all SM particles 
having X = 0. The lightest U{l)x charged state would be absolutely stable (any 
instability would violate X), similarly to the electron, and if also charged under 
U{1)q and/or SU (3)^, an SMP. As above, even if a lighter electrical- and/or colour- 
neutral state with X ^ also exists, the lightest charged one may still be long-lived, 
depending on what decay mechanisms are available. 

A second possibility arises from kinetic mixing between the photon and paraphoton 
which implies that a particle charged only under U{l)x will appear to also have 

a (small) coupling to the photon [951. Experimental constraints on millicharged 
particles [3,98,99] leave a significant parameter space open, including a region of 
interest to accelerator searches, with relatively large charges e ^ 10"'' in units of 
the electron charge and masses m ~ 0.1 — 1000 GeV, with the lower mass bound 
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increasing rapidly above e ~ 10^^. Recently, a 5-dimensional variant of this model 
has been proposed [100] in which the kinetic mixing, and hence the observable 
millicharge, is enhanced as compared to the 4D case. 



In the Standard Model without an additional U{l)x, two additional possibilities 
for millicharged states in principle exist. Firstly, additional SU (3) x SU (2) singlets 
with hypercharge Y — 2edXQ not forbidden, but would be difficult to reconcile with 
grand unified theories [101] and suffer from much tighter experimental constraints 
than in the U{l)x case [99]. A small region of interest to accelerator searches still 
exists at large e = 10^'^ — 10^^ and moderate masses m = 0.1 — 10 GeV. Secondly, 
though massless millicharged neutrinos could be generated by a redefinition of the 
SM hypercharge coupling, this is significantly more difficult in the massive case 
and also implies an unobserved proton-electron charge difference unless e < 10~^^ 
[102]. 



For a conserved non-Abelian gauge symmetry the phenomenon of millicharge is 
excluded and the only possibility is that the lightest state itself carries charge and/or 
colour, in addition to its exotic charge, X. The running coupling constant with Nfx 
exoflavours in SU {N)x is given by the (3 function (see, e.g. [103]): 



_ 1 lliV - 2Nfx 



with /9 < corresponding to a non-confining theory and /3 > to an asymptoti- 
cally free one. For example, QCD is asymptotically free for Nf < 16. Thus there 
are two distinct cases, one f/(l)-like in which the exotic particles are not confined, 
and the second technicolour-like, in which the exotic particles first "hadronise" 
among themselves into composite states of zero total X. Mesonic X-hadrons would 
normally decay rapidly, but the lightest totally antisymmetric SU (N) state (the X- 
proton) could be stable if an analogue of Baryon Number conservation holds in the 
new sector (see, e.g. [104, 105]). A QCD hadronisation would then follow, neutral- 
izing any leftover colour. Some similarity to this picture is found in the hidden- 
valley models of [106], though mainly neutral exotics are considered there. 

Finally, we note that if the new gauge symmetry is broken, then it is hard to see how 
the new states charged under it could be stable, since there should be interactions 
which violate the conservation of the corresponding charge. For the state to be 

long-lived and light, it would have to be an analogue of the muon, which has a 
small mass due to a small coupling with the relevant symmetry -breaking sector 
(the Higgs), but simultaneously a long lifetime due to the gauge boson (the W) 
having a much larger mass. 
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2.4.1 Magnetic Monopoles 



All particles so far observed possess values of magnetic charge consistent with 
zero. The magnetic charge of the electron is constrained by experiment to be < 
^^^^^Qd [107], where go is the elementary (Dirac) magnetic charge. However, in 
spite of a plethora of experimental evidence to the contrary, strong theoretical ar- 
guments continue to motivate searches for magnetic monopoles. 

A potent motivation for searching for monopoles was given by Dirac in 1931 
[108, 109]. Dirac demonstrated that the existence of only one monopole is neces- 
sary to accommodate electric charge quantisation within quantum electrodynamics 
(QED). In a modern form of Dirac's argument [110], the quantisation of the total 
angular momentum in the field of a system of an electric charge q and a monopole 
with magnetic charge g leads to the quantisation condition qg = Here h is 
Planck's constant, c is the speed of light in vacuum, and n is a quantum num- 
ber. Taking g = e as the elementary electric charge, n = 1 sets the theoretical 
minimum magnetic charge which can be possessed by a particle, g^ = n^e. A 
particle with magnetic charge g^ is known as a Dirac monopole. The value of the 
minimum magnetic charge has profound implications for the construction of a the- 
ory of high-energy monopole scattering. The magnetic fine structure constant for a 
Dirac monopole is am ~ rendering perturbative field theory inapplicable. This 
has implications for the reliability of models of monopole production at colliders 
and for the exclusion limits, as described in Sections 4.6 and 7.2, respectively. A 
further consequence is that a Dirac monopole will be largely equivalent to an elec- 
trically charged particle with charge ~ in terms of the electromagnetic force 
it exercises and experiences. This implies that a Dirac monopole will suffer a huge 
electromagnetic energy loss in matter compared with a minimum ionising particle 
(MIP). This is discussed further in Section 5.1.2. 

Although experiments typically search for Dirac monopoles there are a number 
of reasons why alternative values of the minimum charge may be favoured. If the 
elementary electric charge is considered to be held by the down quark then the 
Dirac condition implies that the minimum magnetic charge could be 3gD- How- 
ever, it has been argued that the Dirac condition is not appropriate for a confined 
quark [2, 1 1 1, 1 12] and an observed minimum charge of 3gD may suggest the ex- 
istence of isolated particles with charge |e (free quarks). The existence of isolated 
millicharged particles could imply a minimum charge which is higher still. The 
value of the fundamental charge is also affected if the Dirac argument is applied to 
a particle possessing both electric and magnetic charge. Such a particle is known as 
a dyon [1 13-1 1 SlfJ. It has been speculated that a dyon could exist either as a funda- 
mental particle or as a composite of two particles, one of which possesses electric 
and the other magnetic charge. Schwinger argued that generalising the Dirac condi- 



^ Unless made clear by the context, the term monopole is used to refer to both magnetic 
monopoles and dyons in this paper. 
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tion for a dyon restricts values of n to be even [113-115]. It is therefore important 
that searches are sensitive to as wide a range of magnetic charge as is experimen- 
tally possible. 

While the Dirac argument provides strong motivation for the existence of monopoles 
it gives no prediction of the likely monopole mass. Naive arguments based on the 
classical radius of a Dirac monopole would suggest a mass of the order of a few 
GeV [119]. However, as long as the production cross section is not vanishingly 
small, such a low-mass monopole would long ago have been discovered at collid- 
ers. The available window in mass for Dirac monopoles is largely determined not 
by theoretical arguments but the results of the searches which are described later in 
Section 7.2. 

A further motivation for monopole hunting is provided by their presence in grand 
unified theories (GUTs). 't Hooft and Polyakov showed that monopoles possess- 
ing the Dirac charge, or multiples of it, arise as topological defects of space-time. 
They occur when a simple gauge group is spontaneously broken into an exact U{1) 
subgroup [120, 121]. This would occur, for example, in the phase transition 

SU{5) SU{3) ® SU{2) U{1) SU{3) U{1) (4) 

The spontaneous symmetry-breaking mechanism generates vector bosons X with 
masses, nix- The monopole size Rm can be related to the boson mass via Rm ~ 
rrix^ and the monopole mass is ~ ~ where a is the common gauge 
coupling at unification energy. Here, mx is typically of the order of 10^^ GeV, 
implying a super-heavy monopole of around 10^^ — 10^^ GeV. 

In addition to massive GUT monopoles, it has also been postulated that lighter 
monopoles can be produced through other symmetry-breaking schemes. So-called 
intermediate-mass monopoles (IMMs) with masses between 10'' and 10^"^ GeV 
have been proposed [122, 123]. These could occur in models containing a more 
complicated gauge group than SU (5), such as SO{10). Both GUT monopoles and 
IMMs are beyond the reach of accelerator searches. Such monopoles are sought 
as primordial relics, which could be be bound in matter or found in cosmic rays 
[124-127]. 

To obtain lower-mass gauge monopoles, to which accelerators could be sensitive, 
requires their production via the electroweak symmetry -breaking mechanism. Con- 
trary to earlier work, which asserted that the Weinberg-Salam model could not ad- 
mit monopoles [128, 129], it has been established that monopole solutions are pos- 
sible [130-132]. Furthermore, it has been proposed that monopoles may possess 
masses as low as ~ 1 TeV if the coupling strength of the quartic self-coupling 
of W-bosons was modified, or the group SU{2) U{1) is embedded into a larger 
gauge group [133]. However, as the authors of this work point out, there are theoret- 
ical difficulties with this approach. For example, the modification of the Lagrangian 
may spoil the renormalisability of the Weinberg-Salam model. 
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In higher-dimensional theories such as string theory, space-time is often compact- 
ified on a topologically non-trivial manifold. This naturally leads to a vacuum 
structure suitable for many different kinds of defects, including monopoles. For 
example, in Ref. [134] a magnetic monopole solution is presented in the context 
of heterotic string theory. The mass of this monopole is related to the size of the 
compact space, the string length and coupling. Since these parameters combine to 
determine the low-energy Planck mass and gauge coupling, the mass is constrained 
to be rather large, within a few orders of magnitude of the Planck scale. 

Given the large amount of freedom in higher-dimensional theories with regard to 

the size, dimension and topology of the compact space it is not unrealistic to believe 
that there may be monopole solutions with correspondingly low masses, so that they 
may be discovered at the LHC. 

It should be noted that Dirac monopoles differ from monopoles expected from 
gauge-symmetry breaking. Whilst Dirac monopoles are considered point-like, gauge 
monopoles are expected to have a size i?„j ~ ~ 10^'^^m (for GUTs) and to 
possess a complicated structure of vacuum particles surrounding it. It has been ar- 
gued that, compared to point-like monopoles, the production of gauge monopoles 
from particle collisions will be suppressed by a factor of > 10^° due to form fac- 
tors [135]. 



2.4.2 Q-balls 

Q-balls represent a further possibility of producing topological SMPs. The fastest 
way to understand what they are is to consider a theory of charged scalars with 
SO(2) internal symmetry, in other words a two dimensional internal space of scalar 
fields 01 and 02 with a potential which is only a function of = + (pl- We 
assume that is constant within a (real space) sphere of radius R and zero outside 
that sphere. If (f)i is rotating around the internal S0(2) symmetry with a specific 
angular frequency u, then the conserved charge Q is given by 

/Aqr 
[01^002 - 029o0i] = y , (5) 

whereas the energy E is given by 

E = y i?^ (^..V^ + V) , (6) 

where V is the potential of 0. Using Eq. 5 to replace uj in Eq. 6, it is clear that 
there is a certain radius R at which the energy is minimised, and at this minimum 
E = Q^l^vfcj^. This field configuration or "Q-ball", could decay by emitting the 
basic scalar particle associated with the equations which has unit charge and mass 
777,2 _ Q2y^^02| However the energy per unit charge of the configuration will 
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bQE/Q^ so that if the mass of the scalar m > ^JWJ^, the Q-ball will 

not be able to decay into these scalar particles [136-140]. Many potentials have 
characteristic logarithmic one-loop corrections which allow m > ^2V/(t)'^, not 
least in SUSY. There are examples of Q-balls in both gravity mediated [141] and 
gauge mediated [142] SUSY-breaking models. For generic values of the parame- 
ters in these models, Q-balls must have rather large masses, much higher than the 
electroweak scale, in order to be stable [143, 144]. Although they may be created at 
the end of inflation, it is less likely that they would appear at the LHC. However, if 
present, they could manifest themselves as very highly ionising particles. Coupling 
the scalar field to fermions can make the Q-ball unstable, since decay into those 
fermions would be possible [145]. 



3 Cosmological implications of SMPs at Colliders 

If stable or quasi-stable particles are produced at the LHC, they will also have 
been produced in the early universe. If they are unstable, the lifetime of those par- 
ticles then determines at which cosmological epoch they will decay. Such decays 
generically involve the creation of energetic SM particles, which cause a variety 
of problems that we outline below. In this way one can put constraints upon the 
regions of parameter space where quasi-stable objects might be created. 

If there are relic stable or quasi-stable particles which are produced in the early 
universe then they fall into three broad categories 

• particles which serve as dark-matter candidates. 

• particles which are ruled out since their presence or decay is in one way or an- 
other incompatible with what is observed (see below). 

• particles which exist in such small quantities that they do not serve as dark matter 
candidates and have not yet been detected. 

Of particles which appear in extensions of the SM, very few fulfill the criteria 
necessary to fall into the first category whereas very many fall into the second 
or third. 

"Incompatible with what is observed" can have several meanings, the first being 
due to general relativity. The Friedman-Robertson-Walker solutions of Einstein's 
equations tell us that the rate at which the universe expands depends upon the 
energy density of the contents of the universe. One way a relic particle could be 
incompatible with what is seen is that the presence of that particle would change, 
via gravity, the expansion history of the universe in a way which would be contrary 
to the expansion history derived from observations. The second possible meaning 
is that the particle would decay and its decay products would destroy the observed 
light elements (which are formed in the first few minutes after the Big Bang) or 
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be visible today, for example in the form of gamma rays. A third possible meaning 
is that the particle should have already been detected experimentally by magnetic 
monopole detectors, searches for anomalous isotopes or other experiments of that 
nature. Finally the particle might exist in such an abundance that it gives the cor- 
rect expansion history gravitationally, but its interaction with itself or SM particles 
is such as to rule it out as a dark-matter candidate. 

We will introduce the current state of cosmology before outlining the known prop- 
erties of dark matter. We will then discuss the constraints upon two supersymmetric 
candidates for quasi-stable objects, which could be produced at the LHC, namely 
gluinos and sleptons. Then we will summarise the constraints from cosmology and 
astrophysics upon magnetic monopoles. Finally we will discuss the implications 
for cosmology if a charged massive stable object is produced at the LHC. 

3.1 Cosmology Overview 

Here we will describe briefly how the energy of the universe is divided between 
normal matter, dark matter, and dark energy in a way relevant for the subsequent 
discussion of SMPs at colliders. For more in depth reviews of cosmology see e.g. 
[146]. 

The ratios between the abundances of light elements such as hydrogen, helium, and 
lithium in the universe depend upon the baryon to photon ratio roughly one minute 
after the Big Bang [147]. By comparing their observed abundances with the number 
of photons we observe in the cosmic microwave background (CMB) today, we can 
obtain the number of baryons in the universe. One can also measure the abundance 
of baryons using the CMB alone [148]. 

General relativity provides us with the relationship between the expansion of the 
universe and the density of matter inside it [146]. The density in the form of baryons 
is not enough to explain the observed expansion measured by, for example, the 
Hubble Space Telescope [149]. Furthermore, there appears to be invisible matter in 
galaxies and clusters of galaxies, the presence of which is deduced from the motion 
of gas and stars in those objects [150, 151]. Studies of the clustering of galaxies 
predict that there is 6-8 times as much dark matter as there is baryonic matter [152, 
153]. Observations of type la supemovae tell us that the universe recently started 
to accelerate [154, 155]. 

The conclusions drawn from these different observations are the following - around 
95% of the energy density in the universe is non-baryonic and it couples at most 

weakly to SM particles so that it cannot be seen. Roughly one third of it is dark 
matter while the other two thirds does not dilute quickly as the universe expands, 
i.e. it is dark energy rather than dark matter. At the time of writing, dark energy does 
not appear to be related to particles produced at the LHC in any direct way and we 
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shall not discuss it further. The identity of the dark matter however is obviously 
closely related to the rest of this review. Dark-matter particles may be produced 
at the LHC, so the more we can learn about their behaviour and properties from 
cosmology, the easier it will be to find them. 

If stable massive particles are discovered at the LHC, their couplings and masses 
should not be such that they are over-produced in the early universe. If there is 
too much matter in the universe, then the radiation-dominated regime would finish 
earlier and the universe would be younger by the time the CMB had been red- 
shifted to the temperature that we see today. The universe needs to be at least around 
12 X 10^ years old since observed objects such as globular clusters are at least that 
old. Together with structure formation, this limit provides strong constraints upon 
the over-production of new forms of matter (see e.g. [156]). 

3.2 Constraints on the basic properties and interactions of dark matter. 

Before we look at possible particle candidates for the dark matter in the universe, 
it is interesting to summarise the information we have about its properties. 

3.2.1 Cold dark matter 

By definition, dark matter is at most weakly coupled to Standard Model fields so 
it can only have been produced in the early universe. The simplest mechanism for 
dark-matter production is to have the relics be in thermal equilibrium with the rest 
of the (standard model) plasma at early times. We can use crossing symmetry to 
relate the annihilation rate to the production rate. The criterion for thermal equi- 
librium is that this annihilation rate, which depends upon the coupling and mass 
of the dark-matter candidate, is comparable to the inverse timescale of the expan- 
sion of the universe, which depends only, via gravity, upon the energy density of 
the universe. This means that the universe will not change its size or temperature 
much over the timescale of the particle interaction or decay rate. The dark matter 
observed in the universe today is probably cold, meaning that it was non-relativistic 
at the temperature at which it came out of thermal equilibrium or 'froze' out of the 
plasma. In contrast, the upper bound on the neutrino mass [157] tells us that neu- 
trinos were relativistic until rather late times, at most just starting to slow down at 
matter-radiation equality. Indeed, depending upon the mass hierarchy, some of the 
neutrino species may still be relativistic today. Such fast-moving or hot dark-matter 
candidates will be able to escape all but the largest collapsing structure in the uni- 
verse, so structure formation in a universe with hot dark matter would start with the 
largest structures forming first. 

Observations suggest that the opposite is true, and that the smallest structures form 
first [153] indicating that the dark matter was non-relativistic in the early universe. 
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This provides us with information about the coupling of the particle to the Standard 
Model (annihilation rate) and the mass of the particle: the annihilation rate should 
still be large at temperatures when the particle is not extremely relativistic. 

Strictly speaking however, the fact that dark matter is cold does not mean that it 
had to be non-relativistic at freeze-out; its kinetic energy may have been red-shifted 
away by the expansion of the universe between freeze-out and the epoch at which 
structure formation begins when the universe becomes matter-dominated, a few 
thousand years after the Big Bang. Because of this, one can obtain cold dark matter 
candidates with masses at least as low as the MeV range [158] and probably even 
lower. Most models of cold dark matter under consideration today involve particles 
with masses in the GeV to TeV range, and the combination of this mass range with 
the coupling that will give good relic abundance means that it is possible that they 
could be created at the LHC [159] (one notable exception is the axion [160]). Cold 
dark-matter candidates which interact weakly with Standard Model particles are 
called weakly interacting massive particles or WIMPS. 



3.2.2 Warm dark matter 

Dark-matter haloes in computer simulations have steeply varying density profiles 
in the core, so that typically one would expect that the density of dark matter p oc 
r~'^ where r is the radius from the centre of the halo, and 7 typically has values 
between 1 and 1.5 [161, 162]. However, reconstructions of the distribution of dark 
matter based upon observations of low surface brightness galaxies suggest a much 
less steep density profile. Such galaxies seem to possess an approximately constant 
density core of dark matter [163]. Another discrepancy is that there are many more 
satellite galaxies (smaller galaxies like the Magellanic clouds) produced in A^-body 
simulations of cold dark-matter universes than are observed in nature (for instance, 
in the vicinity of the Milky Way) [164]. There are a number of people who have 
suggested that these discrepancies can be explained if the dark-matter candidate 
is 'warm', in other words half way between cold and hot dark matter [165]. This 
latent thermal energy would smooth out the central density cusps. 

A warm dark-matter candidate produced thermally would usually have a much 
smaller mass than a GeV, so the prospects for dark-matter physics at the LHC 
would change considerably. There is the possibility that weakly interacting neu- 
tralinos are produced non-thermally just before nucleosynthesis [166] so that they 
could be fast moving during structure formation, but producing a relic density of 
neutralinos in this way ignores the fact that one can obtain the correct relic density 
thermally without too much fine-tuning. There are equally a number of people who 
think that these discrepancies between observations and computer simulations may 
be explained using astrophysics rather than changing the properties of cold dark 
matter [167]. 
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3.2.3 Dark-matter self- interaction 



Constraints upon the self-interaction of dark matter can be obtained by looking 
at colliding galaxy clusters. The bullet cluster is an example of two galaxy clus- 
ters colliding, the hot gas forming shock waves at the interface [168]. Gravita- 
tional lensing shows that the two dark-matter populations have moved through each 
other without forming shocks as far as can be observed. This allows one to place 
(mass dependent) constraints upon the self-interaction cross section of dark matter 
of a/M < 2 bam/GeV. This is actually rather a large cross section in terms of 
what we are familiar with in the Standard Model, for comparison agco ~ bam 
(at MeV energies). Some authors have in fact suggested that the dark matter does 
interact with itself, and that this feature solves some of the same problems which 
people try to solve using warm dark matter [169] (see the previous subsection). 
One might therefore think that a colour-singlet particle with a typical QCD cross 
section a ~ ^qcd (for example, perhaps some kind of glueball) might be a viable 
cold dark-matter candidate. However, the constraints on the interaction between 
dark matter and baryonic matter are much stronger (as we shall see in the next 
section). 



3.2.4 Dark-matter coupling to Standard Model particles 

The coupling of the dark matter to SM particles is constrained experimentally in 
two basic ways. The first is by direct experiments, which are looking for momentum 
exchange between the dark matter and sensitive detectors such as CDMS, CRESST, 
DAMA and ZEPLIN [170-173]. Observations of rotation curves of the Milky Way 
suggest that the density of dark matter near the earth is approximately 0.3 GeV 
cm~^ and the velocity corresponding to the gravitational potential of the Milky 
Way is around 200 km s~^. These numbers are based upon the assumption of a 
smooth dark-matter halo and could turn out to be (locally) wrong by more than 
an order of magnitude. Using these numbers, the experimental status is outlined in 
Fig. 3 and the low-energy cross section between a 100 GeV dark-matter particle and 
a nucleon is constrained to be less than ~ 10"''^ cm^. A heavy particle produced at 
the LHC undergoing Standard Model strong or electromagnetic interactions cannot 
therefore be a dark-matter candidate. 

The annihilation of dark-matter particles into Standard Model particles determines 
the rate of interaction of dark matter with the rest of the plasma and hence de- 
termines if the dark matter is in thermal equilibrium with the rest of the universe 
or not [146]. When the dark matter goes out of thermal equilibrium, it essentially 
stops interacting and its number density becomes constant (more precisely the num- 
ber density per co-moving volume - the universe is expanding, so the actual number 
density will decrease due to dilution). 

The analysis of the three-year WMAP data tells us that the density of dark matter is 
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Fig. 3. Constraints upon direct detection of dark matter [170-174]. The red dashed curve is 
from CRESST, purple circles are from Edelweiss 1, green crosses are from ZEPLIN 1 and 
the solid blue curve is from CDMS. The region enclosed by dashes is the claimed detection 
from the DAMA collaboration which does not seem to be confirmed by the other experi- 
ments. The dark green filled region are a set of example SUSY models from Ref. [175]. 



^DAjh'^ = 0.102 ± 0.009 where ^Idm is Pdm/ Pcriu Pent is the density correspond- 
ing to a flat universe [146] and h is the Hubble constant in units of 100 km s^ 
IVIpc^^. In order to obtain this relic density, the thermally averaged cross section 
for dark-matter self-annihilation into Standard IVIodel particles should be [176] 



3 X lO^^W^s"^ 
{av) 



(7) 



so that {av) ~ 3 x 10~^^cm^s~^ is favoured. A cold dark-matter candidate produced 
at the LHC should therefore have this annihilation cross section. This quantity leads 
us to the second method of measuring the coupling of dark matter to Standard 
JVIodel particles, namely through the search for the annihilation or decay products 
of dark matter coming from high-density regions such as the centre of galaxies 
[177]. Since the WIMAP results give us rather good information about {av), the 
uncertainties in this approach lie in the lack of knowledge of the exact density of 
dark matter in dense regions such as the centre of galaxies and in separating the 
signal from dark-matter annihilation from possible background signals. 
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3.2.5 Dark matter at the LHC 



In Section 2 we have learnt about a number of dark-matter candidates. The most 
common SUSY candidates in the literature are neutraUnos. These are WIMPS, 
which are deemed attractive candidates for dark matter since there is considerable 
parameter space in SUSY models where these particles have annihilation cross 
sections rather similar to what is required to obtain the abundance suggested by 
WMAP, as described in Eq. 7 [176, 178]. 

Universal extra dimensions also provide a natural dark-matter candidate if radia- 
tive corrections make the lowest excitation of the hypercharge gauge boson B the 
lightest KK mode [69]. KK modes of other species will decay into the lightest 
mode by emitting zero mode excitations, in other words normal SM particles. The 
lightest mode is stable and if it is aslso weakly interacting it can be a good WIMP 
candidate. There are tight constraints upon the radius and the number of compact 
dimensions if this kind of UED KK dark matter is to provide us with the correct 
abundance, in particular a KK mass around a TeV is favoured, which is encourag- 
ing for LHC studies. If the lightest KK mode turned out to be a particle with EM 
or colour charge rather than the photon then it would be difficult to see how UED 
could remain a good candidate for dark matter. 

There is a large amount of literature on WIMPS (see e.g. [176-179] and references 
therein) and since they are only weakly charged we will not concentrate on them 
here. They do however form candidates for the LSP in SUSY models when the 
NLSP is a charged particle, which then decays into neutralinos. 

3.3 Cosmological constraints on quasi- stable sleptons 

As described in Section 2 if the LSP is the gravitino and the NLSP is the stau, 
i.e. the scalar super-partner of the tau lepton, then the lifetime for the stau decay 
into tau lepton plus gravitino will be relatively large due to the Planck-suppressed 
couplings to the gravitino [180]. 

The gravitino itself is a perfectly good cold dark-matter candidate, although it 
would be impossible to observe through its coupling to the Standard Model fields. 
There is however a potential problem with these dark matter scenarios related to 
nucleosynthesis. 

The epoch of nucleosynthesis in the early universe occurred when the temperature 
was low enough for photodisintegration of nuclei to have ceased but high enough 

for nuclear fusion reactions to occur. Protons and neutrons interchange regularly 
due to beta and inverse beta decay driven by the background thermal bath of neu- 
trinos and electrons. Eventually the plasma becomes cool enough for various light 
elements and isotopes such as deuterium, helium-3 and 4, and lithium-6 and 7, to 
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form. 



Nucleosynthesis gives very specific predictions for the ratio between the light el- 
ements H, He, D, Li etc. that one should observe when one makes observations 
of the abundance of these elements in regions of the universe where they have not 
been processed in stars. 

The strongest constraint upon the quasi-stable stau comes from this observed light 
element abundance. If the decay of the slepton into lepton plus gravitino takes place 
during or after nucleosynthesis then energetic particles will be injected into the 
plasma, leading to various effects such as the dissociation of these light elements. 

In particular, the injection of energetic protons into the plasma due to the decay of 
unstable relic particles will lead to the reaction p -|-^He — >^He -I- D which increases 
the abundance of deuterium [181, 182]. This deuterium can also subsequently in- 
teract with helium to increase the abundance of ^Li. 

The ^He/D ratio is a strong function of the decay of the slepton into gravitino 
[183, 184] and the observed value of ^He/D restricts the stau mass ruf to be greater 
than about Sm^ for rUf ~ 500 GeV or 20mg for ruf ~ 100 GeV. This cuts into the 
interesting region for collider searches. For instance, one can in principle recon- 
struct the gravitino mass by simply measuring the mass of the stau and the energy 
of the emitted tau [180] 

m| = m? - 2mfEr (8) 

but m^^ can only be measured in this way if it is not too small. Since the constraints 
from light element abundance restrict the gravitino mass for a given stau mass, 
much of the interesting parameter space is ruled out. One way of avoiding this 
is for there to be entropy production at some temperature below the freeze-out of 
staus from the plasma and the era of nucleosynthesis. Such entropy production will 
increase the number of relativistic particles in the plasma and effectively dilute the 
sleptons, so that when they do decay, their effect upon the light element abundances 
is diminished. We will mention this more in Section 3.5. 

One can also consider situations where the superpartner of the axion, the axino 
is the LSP. This gives rise to analagous dark matter physics, although different 
constraints from nucleosynthesis [185]. 



3.4 Stable and quasi-stable gluino 

In this section we present the cosmological constraints upon the stable and quasi- 
stable gluino. 
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Fig. 4. Constraints upon the gluino mass and SUSY breaking scale parameters in models of 
split supersymmetry. Shaded regions are ruled out by different cosmological observations. 
The tightest constraints come from the diffuse gamma ray background for rrig < 300 GeV 
and from Big Bang nucleosynthesis (BBN) for nig > 300 GeV. The plot is taken from 
Ref. [182]. 

3.4.1 Stable gluino 

The case of the stable gluino is presented in detail in Ref. [26]. If the gluino is 
the LSP then the relic abundance is calculated by obtaining the thermally aver- 
aged annihilation rate and comparing it to the expansion of the universe, giving the 
approximate result 

2 



(9) 



VlOTeV/ 

The stable gluino could not be the dark matter as its coupling with Standard Model 
fields is greatly in excess of cross sections that are already ruled out (see fig. 3). 
We therefore require Qgh"^ <^ 0.1, giving a maximum mass for the gluino of a few 
TeV. These gluinos will be captured by nuclei, leading to the possible detection of 
(highly) unusual isotopes of nuclei [186]. Non-observance of these isotopes restrict 
values of > 100 GeV although non-perturbative effects could change the values 
predicted by equation (9). 



3.4.2 Quasi-stable gluino 

The constraints upon the quasi-stable gluino come from the same places as the 
constraints upon the slepton described in the previous section. We have seen in 
some versions of split supersymmetry [46] that gluinos can be produced which 
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then decay outside the detector. The Ufetime of these gluinos is given by [48] 



where ms is the scale of supersymmetry breaking, which is typically very large in 
these models, » 10^° GeV. The gluino lifetime as function of ms is displayed in 
Fig. 2. As for the stable gluino, constraints from the non-observance of anomalous 
isotopes become stronger as one raises the energy of ms, since the anomalous 
isotopes will be increasingly stable for cosmological time-scales. Fig. 4 shows how 
different cosmological constraints rule out values of ms which depend upon mg, 
for example, ms must be less than 10^^ GeV for a gluino mass mg ~ 200 GeV or 
less than 10^^ GeV for m~g ~ 10 TeV [182]. 

As one lowers the value of ms, the lifetime of the gluino falls, and the annihilation 
products become progressively more troublesome for cosmology. As one would 
expect, as the lifetime drops, the time at which the strongest constraints emerge 
relates to progressively earlier epochs in cosmology. For example, at values of ms, 
lower than those ruled out by searches for heavy isotopes, the parameter space is 
ruled out through constraints upon the diffuse gamma-ray background today made 
by the EGRET gamma ray observatory [187, 188]. 

Values of ms greater than about 10^^ GeV are ruled out since the decay of such 
gluinos would lead to distortions in the thermalisation of the CMB before recom- 
bination. The CMB is the radiation remaining from the time when the temperature 
of the universe decreased such that the mean free path for photons approaches in- 
finity. At this epoch, a few hundred thousand years after the Big Bang, the universe 
becomes transparent. The CMB is the most perfect blackbody ever observed, so 
distortion of the spectrum by the injection of high-energy photons into the plasma 
leads to strong constraints [189]. 

Finally gluinos with very short lifetimes will, like the unstable stau, cause problems 
for light-element abundance. The decaying gluinos will photo-dissociate the light 
elements left over after nucleosynthesis although for mg < 300 GeV the gluino 
lifetime is long enough to avoid the nucleosynthesis constraint [182]. 

Very roughly, the combined constraints upon the gluino are as follows - if the gluino 
mass mg < 300 GeV then the SUSY scale ms < 10^^ GeV whereas if mg > 300 
GeV then ms < 10^° GeV. 




(10) 
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Fig. 5. Amount of entropy which needs to be produced in order for the unstable stau to be 
consistent with nucleosynthesis constraints, see Eq. 11. Plot taken from [183]. 

3.5 Consequences of the observation of a charged stable object at the LHC 
3.5.1 Electrically charged particles 

If a charged massive stable object is discovered at the LHC then it may be very 
interesting for cosmology. If the object were truly stable, then it would be difficult 
to understand why it is not observed today. We have already discussed that such a 
charged object could not be the dark matter. On the other hand, if the object were 
found to possess a lifetime such that one would expect its decay to be problematic 
for nucleosynthesis or the CMB, or that it would overproduce gamma rays, again it 
would be difficult to understand. 

If such a situation were to arise, it would be necessary to invoke a mechanism, 
which would dilute the density of the relic. Let us call our problematic charged 
stable particle X. If there were another species of particle Y , which decayed in 
between the time at which the stable relic X froze out of the plasma and the time at 
which it {X) causes problems (by decaying or by not being observed), then it might 
be possible to dilute the particle. The decay of the particle Y may produce more 
light particles such as photons, thereby reducing the number density per photon 
of the dangerous relic X. Another way of looking at this is to say that the decay 
of particle Y into light particles temporarily stops the temperature in the radiation 
bath from dropping as the universe expands. At the same time, the number density 
of particle X would drop as it is no longer in equilibrium with the radiation bath 
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but has already frozen out. In each particular scenario, the dangerous relic particles 
would have to be diluted by a factor A in order for their decay not to be dangerous 
for light element abundance. The dilution of the dangerous X particle (in this case 
the stau) is written 

IT-Xbefore _ 1 f^Xafter ^-^y^ 

s As 

where nx is the number density of the X particle and s is the entropy density, 
which is roughly equal to the number of relativistic particles in the plasma per 
unit volume. As an example. Fig. 5 shows the amount of entropy, written as the 
required value of A, which should be produced after freeze-out in order for the 
unstable stau described in Section 3.3 to be compatible with constraints from light 
element abundance [183]. 

Another effectively equivalent way of diluting the density of a stable relic is via 
inflation and reheating. In principle there could be a period of inflation at an energy 
scale lower than that of the LHC. However, if this period of inflation were driven 
by an inflaton it would have to be rather strongly coupled to the Standard Model 
particles in order to achieve thermalisation before the beginning of nucleosynthesis 
(see however [190-192]). 



3.5.2 Magnetic Monopoles 

Theories which contain the possibility of spontaneous breaking of some gauge sym- 
metry to a true vacuum with a non-trivial topology contain solitonic objects, which 
are topologically stable. These objects are simply classical field configurations of 
the equations of motion. The simplest example of this is when there are two discrete 
degenerate vacua; this is the familiar case of a massive real scalar field with self- 
interaction and a negative mass-squared term. In separate, causally disconnected 
regions in the early universe, the field will take a random decision as to which 
vacuum it will inhabit when the temperature drops below the critical temperature 
corresponding to the phase transition. Later, as the horizon grows and different 
parts of the universe come into causal contact with each other, regions are found 
where the field is forced to interpolate (smoothly due to the equations of motion) 
between one vacuum and another vacuum. This region corresponds to a domain 
wall, a wall trapped in the symmetric state of the early universe. For different vac- 
uum structures, different classes of objects can be produced such as strings for a 
U(l) vacuum. It is in this way that monopoles are formed cosmologically, when the 
vacuum state is a continuum of degenerate vacua with a spherical topology in field 
space. For more details on different topological objects, see [193, 194]. 

The idea that there should be on average one monopole created per causally dis- 
connected region in the early universe was first intoduced by Kibble [193]. Any 
theory of particle physics where there is too much non-relativistic matter created 
in the early universe gives problems for cosmology since no monopoles have been 
detected. One of the initial motivations for the theory of inflation was to provide a 
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dilution of the density of monopoles created in the early universe [195]. 

A very tight constraint upon the density of monopoles comes from the observa- 
tion by Parker that the presence of magnetic monopoles would limit the ability of 
astronomical objects to build up magnetic fields [196]. The constraint for interstel- 
lar space is that there should be less than one monopole per 10^^ nucleons today. 
This means any symmetry-breaking scale larger than 10^° GeV would lead to an 
unacceptably large number of magnetic monopoles if inflation does not occur. 



4 Modelling the production of SMPs at colliders 

In this section the techniques used to model the production of heavy particles at 
colliders are discussed. The aim of this section is not to provide an overdetailed de- 
scription of the phenomenology of proposed, undiscovered particles, but instead to 
show how techniques and models developed for SM particles can profitably be used 
to predict the gross features of SMP production at colliders. We focus mainly on 
SMPs which do not possess magnetic charge in view of the inapplicability of per- 
turbative field theory to monopole processes. A description of monopole production 
mechanisms which have been considered in collider searches, and the approxima- 
tions which have been used to calculate their rates is given at the end in Section 
4.6. 



4.1 Production rates 

Exotic stable or long-lived new particles are usually thought to be pair-produced at 
a collider: 

ah ^ X,Xd, (12) 

where a and h are normal constituents of the incoming beams, and Xc and X^ be- 
long to the same new theory. While pair production may occur at a large rate, the 
decay of particles X^. and X^ may be suppressed by the existence of a new (almost) 
preserved quantum number, possibly resulting in long-lived particles which could 
interact in the detector as an SMP. Besides direct production, SMPs could be pro- 
duced via the decay of a heavier particle, for example if a state Xi would decay 
rapidly into a lighter state, which in turn cannot (rapidly) decay any further. 

If the production of a particle X is allowed singly 

ah X, (13) 

then the decay X ^ ab may also occur. The production cross section is directly 
proportional to the decay width, a{ah — > X) oc r{X — > ah), where the constant of 
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proportionality involves mass, spin and colour factors, as well as the ab partonic 
flux, but no model- specific details. A singly-produced particle which is sufficiently 
long-lived so as to penetrate through a detector, has, therefore, a production rate 
which will be negligible, as can be seen from the following example. Consider a 
colour-singlet scalar state with mass 200 GeV, coupling only to gluons. A lifetime 
of cr = 1 mm will translate into a small production rate of 0.0002 events per year, 
assuming full LHC luminosity, 100 fb^. The larger the lifetime, the smaller is the 
production rate. If the state is a colour octet, there is a factor eight enhancement, or 
if the particle has spin one, another factor three must be added. However the produc- 
tion rates will still be negligible. A potential exception could be off-shell decays, 
such as X — > W'^W~ for mx <^ 2mw, where F may be strongly kinematically 
suppressed, while a would be less affected since the incoming W's, emitted off 
the beam constituent quarks, are spacelike. In practice, taking into account current 
experimental limits discussed in Section 7 as well as the impact of loop-induced 
decays, such prospects are excluded. 

Given the Lagrangian of the theory, and the values of relevant couplings and masses, 
the parton-level differential cross section for da{ab — > XcXd) is readily obtained. 
Convolution with the a and b parton densities gives all relevant differential distri- 
butions, and integration over phase space provides the total cross section. The pair 
production cross section of coloured particles is of 0{al), while that of colour- 
singlet ones is of 0{al,^), where a,.^ represents electroweak couplings. Thus the 
production rate for colour-singlet particles would be a factor {aew/ «s) ^ times smaller. 
Additionally, the latter cannot couple directly to the gluonic content of the beams, 
so typical production rates for colour-singlet particles are about 2 to 4 orders of 
magnitude lower in rate. To illustrate this. Fig. 6 shows Pythia predictions for the 
pair-production cross section at the LHC of exotic fourth-generation quarks with 
charge ±|e and fourth-generation leptons with charge ±e, as a function of the 
mass of exotic particles. The lepton cross section is clearly dwarfed by that of the 
quarks. However, this argument only concerns the direct production mechanisms. 
Weakly interacting particles can be produced at a large rate in the cascade decays 
of heavier, strongly interacting particles. 

Next-to-leading-order (NLO) calculations of event rates have been performed for 
many exotic scenarios, such as SUSY [197-199], but not for all. Experience shows 
that K = crNLo/<7LO typically range between 1 and 2 [200]. In this article K = 
1 will be assumed throughout, which will provide conservative estimates of the 
experimental feasibility of discovery. 



4.2 Event topologies 

Higher-order perturbative calculations involve emissions of further partons. This 
can be understood as ab — > XcXd generalising to ab — > XcXae and ab — > 
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Fig. 6. Pythia predictions of the pair-production cross section for fourth-generation 
quarks of charge ±|e, and leptons of charge ite at the LHC. 

XcXd ef, etc., where e and / will mainly be gluons for incoming hadron beams, 
but also could be quarks, leptons or photons. Fig. 7. Such emissions modify event 
shapes and thereby the experimental signatures. 

These processes can be described either by higher-order matrix elements or by 
parton showers applied to lower-order matrix elements. There are relative advan- 
tages and limitations for both approaches. The former are more accurate for a few 
well- separated emissions, and of course contain all the model-specific details, but 
they diverge in the soft/coUinear limits and are not meaningful there. The latter 
are based on model-independent approximations that work especially well in the 
soft/coUinear regions, in which dampening by Sudakov factors and resummation of 
multiple emissions ensures a physically meaningful behaviour. In practice, showers 
have turned out to be reasonable approximations up to the scale of the lowest-order 
process, i.e. typically for the additional-jet transverse momenta below the mass 
scale of the produced X particles [201]. For exploratory studies the shower picture 
alone is therefore sufficient, while, for precision studies, matrix-element informa- 
tion also has to be supplied. Various methods to combine the two approaches, so 
that the matrix-element behaviour is reproduced at large separation scales and the 
parton-shower one at small scales, are under active study [202-206]. Following a 
discovery, both higher-order matrix-element calculations and corresponding match- 
ings would have to be carried out. 

In the shower approach, one distinguishes between initial-state radiation (ISR) from 
the incoming a and b partons and final-state radiation (FSR) from the outgoing 
Xc and Xd, with interference effects neglected. In the current case, since the Xi 
particles are heavy, radiation from them is strongly suppressed [207] in comparison 
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Fig. 7. Top: examples of production of squarks (top left), gluinos (top center) and heavy 
leptons (top right) in hadron-hadron colliders. Bottom: examples of production of squarks 
(bottom left), gluinos (bottom center) and heavy leptons (bottom right) in e"'"e~ -colliders. 
Examples of ISR and FSR are displayed. 

to light particles. Thus ISR dominates. The main experimental consequences of ISR 
is that the XcX^ pair is produced in association with a number of further jets, and 
that Xc and X^, do not have opposite and compensating transverse momenta, as 
they do in the LO picture. This p±_ imbalance is in itself a revealing observable, 
especially when the Xj are not directly detectable. On the other hand, the activity 
of the additional jets, and of the underlying event, may be a nuisance for some 
studies. Jets from FSR could be important, as will be discussed in the following 
section. 

The kinematics of heavy particles produced in collisions can be studied with the 
help of several programs. The two general-purpose generators Pythia [208] and 
Herwig [209, 210] are traditionally used to study physics within and beyond the 
Standard Model. They are mainly based on leading-order matrix elements, but in 
a few cases also NLO matching is available. Their strength is that they provide 
complete event topologies, such as they could be observed in a detector, includ- 
ing descriptions of parton showers, underlying events, and hadronisation. Many 
of these aspects are handled in different fashions in the two programs, e.g. with 
respect to the choice of shower evolution variables or hadronisation schemes, but 
they tend to give similar results. They have both been supplemented with routines 
to handle hadronisation into SMPs, as will be described in the next subsection. 
Other general-purpose generators, but currently without hadronisation into SMPs, 
are Isajet [44] and Sherpa [211]. 

Many other programs provide matrix elements to LO or NLO. A comprehensive 
survey is given in the "Les Houches Guidebook" [212], while dedicated BSM tools 
descriptions and an online repository can be found in [15, 16]. Examples of pro- 
grams that can calculate arbitrary LO processes, once the Feyimian rules have been 
encoded, are MadGraph/SMadGraph [213-215] and CompHep/CalcHep 
[216,217]. Examples of programs that contain NLO matrix elements are Prospino 
[218] for SUSY particle production and SDecay [219] for SUSY decays. In UED, 
production cross sections for gluon and quark KK excitations at hadron colliders 
have been calculated [66] and are available as an extension [75] of the Pythia 
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generator. These programs do not provide complete events, however, and therefore 
have to be interfaced with general-purpose generators to be fully useful. This pro- 
cess is rather well automated: the original Les Houches Accord (LHA) [220,221] 
provides a standard format to transfer simple parton-level configurations, includ- 
ing information on flavours, colours and momenta, while the SUSY Les Houches 
Accord (SLHA) [17] allows standard exchange of SUSY parameters, couplings, 
masses, and branching ratios. 



4.3 Hadronisation 



For colourless new particles the story ends here. Once produced, the Xj's will sail 
out towards the detector. A coloured particle, however, cannot escape the confine- 
ment forces. Therefore it will pass through a hadronisation stage, during which it 
picks up light-quark or gluon degrees of freedom to form a colour-singlet "hadron". 
For a colour-triplet, denoted C3 to highlight its coloured nature, this may either be 
a "meson" C-iq or a "baryon" Csqiq2. For a colour-octet, denoted Cg, the alterna- 
tives are a "meson" Csqiq2, a "baryon" C^qiq2q^ (or corresponding antibaryon), or 
a "glueball" C^g. Collectively, such states will be referred to as i?-hadrons, irre- 
spective of the physics scenario which gives rise to them. The name is borrowed 
from SUSY, where it refers to the nontrivial i?-parity possessed by such hadrons. 
Specific states will be denoted by an R with the flavour content as lower indices and 
the charge as upper one[ll e.g. for a gud state. In order to avoid too unwieldy 
a notation, spin information is suppressed. 

The hadronisation process does not appreciably slow down the _R-hadron relative 
to the original C3 g, as can be seen from the following. Consider a colour triplet C3 
with energy Ec^ and longitudinal momentum pycg (here longitudinal is defined in 
the direction of the C3), which hadronises into an i?-hadron with energy En and 
longitudinal momentum py^^, plus a set of "normal" hadrons that (approximately) 
take the remaining energy-momentum. These hadrons are produced in the colour 
field pulled out behind the C3, and therefore ought to be limited to have smaller 
velocities than that of the endpoint C3 itself or, more precisely, to have smaller 
rapidities y defined with respect to the C3 direction. Now recall that E+p\\ = m^e^, 
where m± is the transverse mass of a particle. So if a normal hadron h, with a 
{m±h) ~ 1 GeV, could at most reach the same rapidity as the C3, then it follows 
from the above that the ratio between the E + value of the normal hadron and 
that of the C3 could at most be {m±h) / mc^ ^ 1 GeV/mcg. Finally, assume that 
the normal hadrons are produced behind the C3, located at some rapidity yo, with a 
typical rapidity separation of (Ay) 0.7 In 2, as in ordinary jets, i.e. at yo — In 2, 

?/o — 2 In 2, I/O — 3 In 2, The above upper bound for E +p^^ of a single hadron then 

is replaced by an average for the summed effect of the normal hadrons produced in 



This is the convention adopted in the 2006 Review of Particle Physics [1]. 
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association with the C?;. m±_ cxp(7/o)(exp(— In 2) + exp(— 2 In 2) + exp(— 3 In 2) + 
■ ■ ■ ) = cxp(yo), where each term expresses the E+p^^ value of the next hadron. 
The i?-hadron will retain all the C3 energy not taken away by these normal hadrons. 
One therefore arrives at the expectation [222] that the fragmentation parameter z 
typically is 

(.)h|^>m|^)«i-^-i-^. (14) 

Ec3+P\\C3 mcs 



That is, if C3 has a large mass, (z) is close to unity, and all other normal hadrons in 
the jet take a very small amount of energy, ~ 1 GeV • 7^3, where 7^3 = Ec^/mcs 
is not large for a massive state C3. For a colour octet the energy loss from the 
i?-hadron to the rest of the jet would be about twice as big as for a colour triplet 
(Ca/Cf = 9/4), but otherwise the argument for colour triplets applies in the same 
way. 

The above z value is not accessible experimentally since, in a busy hadronic en- 
vironment, it is not possible to know which particles come directly from the C3 § 
hadronisation. Instead, with an i?-hadron being part of a jet, the natural observ- 
able is the fraction z* of the jet energy, Ejet, which is carried by the i?-hadron: 
z* — Eji/Ejet- Fig. 8 shows the expected distribution of z* for gluino i?-hadrons of 
mass 300 GeV produced at the LHC. The jets were reconstructed with the CDF Run 
I cone algorithm with R = 0.7 [223]. Predictions are shown using the string and 
cluster models of Pythia and Herwig, respectively. The details of these models 
are given later in this section. Also shown are Pythia predictions for the distri- 
bution of z* of the leading particles contained in jets produced in Standard Model 
2-to-2 QCD processes (labelled SM-jets). 

Note that the i?-hadron (z*) is significantly smaller than the prediction of Eq. 14, 
(z) ^ 1 — 2/300 0.993, including a factor of two for the colour-octet nature of 
gluinos, but not taking into account the possible loss of particles outside the cone. 
Instead several other physics components contribute to the final curve, by deposit- 
ing further particles inside the jet cone. 

• Final-state radiation from the gluino, although small, is still larger than the non- 
perturbative effect discussed above. Since coUinear FSR is strongly suppressed 
[207], the additional FSR jets will not be correlated with the gluino direction, and 
so may or may not end up inside the R-jet cone. At hadron colliders the FSR jets 
tend to drown among the more numerous ISR jets. They could stand out at a lepton 
collider, at which there is no QCD ISR. 

• Initial-state radiation from the incoming quarks and gluons is not suppressed by 
any mass effects. The ISR jets may appear anywhere in phase space, by chance also 

inside the R-jet cone, and may then deposit a significant amount of extra energy. 

• The underlying events ("beam remnants" and "multiple interactions") tend to de- 
posit particles inside any jet cone, in proportion to the size of the cone. Normally 
these particles would be rather soft, but upwards fluctuations can occur. 
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Fig. 8. The predicted distribution of the fragmentation variable z* = Er/E-^^x for gluino 
i?-hadrons of mass 300 GeV produced at the LHC. The expectations of the string and 
cluster hadronisation schemes, implemented in Pythia and Herwig, respectively, are 
shown. Also shown is the prediction from PYTHIA of the z* distribution of leading particles 
within SM jets. 

For the actual hadronisation of the coloured heavy object there are two main models 

available. In PYTHIA the Lund string fragmentation model [224] is used, wherein 
an assumed linear confinement potential is approximated by a string with a con- 
stant tension of k pa 1 GeV/fm. A colour triplet Cs (or C3) is at the endpoint of 
such a string. When the C3 moves away from its production vertex it pulls out the 
string behind itself, to which it loses energy. This string may then break by the 
production of a light quark-antiquark or diquark-antidiquark pair, where the q or 
qiq2 is in a colour-antitriplet state which can combine with the C3 to form a sin- 
glet. Further breaks of the string may occur, which causes the formation of a jet 
of hadrons. The composition of the different light flavours is assumed universal 
and thus constrained e.g. by LEP data [225, 226] (u : d : s ^ 1 : 1 : 0.3, with 
diquarks further suppressed). A mixture of different charge states is finally pro- 
duced. The energy-momentum fraction z retained by the i?-hadron is described by 
an appropriate fragmentation function, with parameters that fit e.g. B meson spec- 
tra [227,228]. These functions have a mass dependence consistent with Eq. 14 for 
the extrapolation from b to C3 hadronisation. 

A colour octet Cs in the Lund string model is viewed as the incoherent sum of a 
colour and an anticolour charge (the planar or A^^ 00 approximation [229]), 
such that an octet Cg is attached to two string pieces instead of one only. Each of 
these pieces is allowed to break as above. One gives a quark or an antidiquark, 
the other an antiquark or diquark, and these combine to form an i?-meson or R- 
baryon. Diquark-antidiquark "hadrons" are rejected. In addition, a new mechanism 
is introduced: a gg pair may be produced, such that an i?-glueball is formed and 
the leftover gluon attaches to the two string pieces. Its relative importance is not 



38 



known. The equivalent mechanism for nomial gluons would lead to the production 
of glueballs. There are a few indications that this may occur, [230, 231], but it 
cannot be at a significant rate. Based on the absence of real evidence for normal 
glueballs [232], the default value for the it!-glueball fraction in Pythia is 10%, and 
can be changed. Clearly a scenario in which glueballs are predominantly formed in 
the hadronisation step will have a large impact on experimental searches which 
rely on the reconstruction of tracks. However, as is outlined in Section 4.4, the 
behaviour of an i?-glueball when it interacts in material is expected to be similar 
to that of a neutral i?-meson and it may therefore convert into a charged ii!-hadron 
which would leave behind a track. 

Predictions from Pythia of the fractions of different i?-hadron species formed in 
the hadronisation of a gluino and a stop are shown in Tab. 3. The predicted species 
of an it!-hadron arising from an antistop are almost exactly the charge conjugates 
of a stop i?-hadron. The key observation is that, when neglecting the a priori frac- 
tion of i?- glueballs of 10%, roughly 50% of the produced states are charged. The 
gluino /?-baryon numbers are somewhat lower than expected. This is an artifact of 
approximations used in hadronisation model. However, these should anyway con- 
tribute less than 10% of the i?-hadrons produced. It is also interesting to note that, 
as discussed in Section 5.2, although i?-hadrons will largely start out as mesons, 
nuclear scattering in calorimeters will ensure that they mostly end up as baryons as 
they leave the detector. 

Herwig is based on cluster fragmentation. In this picture all gluons are nonper- 
turbatively split into quark-antiquark pairs at the end of the perturbative cascade. 
Partons from adjacent such breaks, and from original (anti)triplets, are then grouped 
into colour- singlet clusters. Large-mass clusters are broken into smaller ones, along 
the "string" direction. These clusters then decay to two hadrons, using phase-space 
weights to pick between allowed flavours. A colour octet Cg is combined with a 
gi and a from two adjacent nonperturbative splittings to form a cluster [233]. 
The cluster decay can be either of Csqiq2 — ^ Csq-iq2 + Csqi% + or 

Csg + qiq2- The two former possibilities are handled as usual based on phase 
space weights, while the latter i?-glueball possibility is added with a free normal- 
isation. By default it is set to zero. The resulting gluino i?-hadron composition is 
shown in Tab. 3. 

Thus, although the technical details in Pythia and Herwig are quite different 
the basic principles are similar. In both models, i?-hadron production occurs pre- 
dominantly by the C^^g picking up u and d quarks, while s is more rare, and c or 6 
can only be produced in the shower, at an even lower rate. Apart from the unknown 
fraction of R-glueballs, the u-d charge difference ensures that approximately half of 
the produced _R-hadrons will be charged. The two generators should therefore pro- 
vide similar phenomenology for the flavours and momenta of produced /?-hadron 
events. Any differences would reflect true uncertainties in our current understand- 
ing. 
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Table 3 

Predictions from PYTHIA of the fractions of different species of i?-hadrons following the 
hadronisation of a gluino (left) and a stop (right) of mass 500 GeV produced at the LHC. 
The Herwig gluino predictions [233] are for a 2000 GeV mass, but almost identical for 
50 GeV. 

4.4 R-hadron properties 



The mass splittings of i?-hadrons are of critical importance when designing a search 
strategy. If one state would be significantly lighter than another, one would expect 
this state to be dominantly present in the detector. If the lightest state would be 
neutral, a completely different (and far more experimentally challenging) signature 
is expected than if the lightest state would be charged. 

The masses of the produced i?-hadrons are best understood from the mass formula 
for the lowest-level (i.e. no radial or orbital excitation) hadrons [234, 235] 

E, (Fj • Fj) (Si • Sj) 
rrii-k}^ — (15) 
. ^ rui ruj 

where rrii are the constituent masses, Fj are the colour SU(3) matrices, Sj the spin 
SU(2) ones, and k a parameter related to the wave function at the origin. The C3.8 
is so heavy that it provides an almost static colour field in the rest frame of the 
i?-hadron, and therefore its spin, if any, is decoupled (cf. "Heavy Quark Effective 
Theory" [236]). The heavy (73 § has thus a strongly localised wave function, while 
the light degrees of freedom are spread over normal hadronic distance scales. The 
relative localisation of the wave functions largely accounts for the expected en- 
ergy loss and scattering behaviour of i?-hadrons in matter, as discussed in detail 
in Section 5.2. It also ensures that the mass splittings, given by the second term in 
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Eq. 15, are determined by the light degrees of freedom. It must be noted that the 
"experimentally observable" constituent masses used below are different from the 
renormalisation-scheme-dependent running masses found in the Lagrangian of a 
theory although it is possible to provide an approximate translation for any specific 
case. The mass splittings discussed below agree with those obtained earlier with a 
similar bag model approach [237,238], and with lattice calculations [239]. 

For a colour-triplet C3 the hadron mass is easily obtained from standard quark 
and diquark constituent masses. For C3 mesons (CsQ, C^q), no significant mass 
splitting is expected to occur [235], in analogy with the absence of mass splittings 
of B-hadrons. For C3 baryons such as C'^qq states (and, by symmetry for the Cggg 
baryons) the mass spectrum would be (all units in GeV) 



(-1 X 

Mc.qq ^ Mc, + 0.3 + 0.3 - 0.026 x \l ^^J s,. 



0.3 X 0.3 



;Mc3 +0.46 



Mc3,, ^ Mc, + 0.3 + 0.3 - 0.026 x ^^^^^^^ 



^ + 0.65 
where Sqg denotes the total spin of the qq system. 

In the case of a Cg state, i?- mesons, i?-baryons and i?-glueballs could arise from 
hadronisation. The i?-mesons Csqiq2 give a colour factor Fi ■ F2 = 1/6 for the 
light quarks, to be compared with —4/3 in a normal ^1^2 meson. For a Csqq state, 
the mass spectrum is thus given by: 



Mcsgq^Mc, + 0.3 + 0.3 - 0.043 x ^^^^ 



■ + 0.66 



ix+i) 



Mcsqg ~ Mcs + 0.3 + 0.3 - 0.043 x f^^^ s^q = 1 

Mcg + 0.58 

The p-TT mass difference therefore flips sign and is considerably reduced in size 
for the corresponding it!-mesons. Further, a gluon constituent mass is about twice a 

light-quark [240,241], so the C^g state is almost mass degenerate with the C^uu and 
Csdd ones, and will have similar properties. Assuming that the constituent mass of 
a gluon is approximately 700 MeV [240,241], the mass of a C^g state is then 

Mc,g ~ Mc, + 0.7. (16) 

The expressions for the masses of the i?-baryons (Cgqqq) are somewhat more cum- 
bersome, since there are more colour and spin combinations possible, but do not 
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offer any peculiar features. Splittings are expected to be less than the order of a 100 
MeV [235], and the mass of the Cgqqq states is thus approximately 

Mc,,,q = Mcs + 0.9. (17) 

Furthermore, the orbital and radial excited states would follow patterns not too dis- 
similar from normal hadronic spectra. However, before a discovery of it!-hadrons 
it is not useful to attempt to calculate the full hadron spectrum including excited 
states. In case such an excited state is produced it would probably decay electro- 
magnetically or strongly at such a short time scales that no secondary vertex is re- 
solved. Owing to the small mass splittings and the small boost factor of a massive 
i?-hadron, the additional photon or pion produced will also be of low momentum, 
and so drown in the general hadronic environment. For the current studies it is suf- 
ficient to simulate the production of the lowest-lying state of each allowed flavour 
combination. 

Although this paper is only concerned with stable i?-hadrons, it is worthwhile to 
briefly consider the possibility of 7?-hadron decays. If unstable over nanosecond 
time scales, i?-hadrons could decay inside the detector. This would take place as an 
almost free decay of the Cs,8, with the rest of the hadron acting as spectators, and 
could therefore be described by the standard perturbative picture. Showering and 
hadronisation can be added, the latter also involving the spectators. These decays 
may well violate baryon or lepton number, via such processes as ^ — > uds, and give 
rise to unusual hadronisation scenarios [242, 243]. 

A further property of it!-hadrons about which little is known is their oscillatory 
behaviour. Neutral i?-mesons will be able to convert into their anti-particles without 
violating any known conservation laws. Since the oscillation length depends on 
couplings and masses of particles which have yet to be discovered (if ever) then 
either minimal or maximal oscillations are conceivable [244, 245]. As discussed 
in Section 5.2 this may have implications for the ability of experiments to both 
discover and quantify a heavy squark which may be produced at a collider. 

4.5 Sources of uncertainties 

The understanding of the production of non-magnetically charged SMPs is not 
complete, but should be fully adequate for the purposes of a search-and-discovery 
mission. Further sophistication would be added once the first signals for new physics 
would point the way to a more specific scenario, rather than the generic ones consid- 
ered here. In any search relying on QCD models, the following areas of uncertainty 
may be relevant: 

• Production rates: are very much model-dependent, and contain some uncertain- 
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ties from uncalculated higher-order contributions. These can be parameterised in 
the form of a fT-factor. 

• Event topologies: contributions from further parton emissions are treated by the 
parton-shower approximation, which should be sufficient for the bulk of the cross 
section, but not necessarily for topologies with well-separated further jets. 

• Hadronisation, momenta: the fragmentation function is guaranteed to give most 
of the momentum to the /?-hadron, so uncertainties are restricted to the minijet of 
a few further normal hadrons produced in the same general direction. The typical 
minijet energy could easily be uncertain by a factor of two. 

• Hadronisation, flavours: some uncertainty, especially in the glueball and baryon 
sectors. The rate of the former is a completely free parameter. It should not affect 
the bulk of the production, mesonic states with u and d quarks only, which also 
will guarantee a rather even mixture of charged and uncharged states. 

• i?-hadron properties: the general pattern of states and masses appears to be well 
understood, especially in the dominant meson sector. Should, against all expec- 
tations, decays like Csud C^g + 7r+ be kinematically allowed, they would 
have consequences for the charged/neutral ratio. 

4.6 Production mechanisms of for magnetic monopoles 

The interactions of magnetic monopoles at high energies are difficult to treat the- 
oretically since perturbation theory is inapplicable due to the size of the elec- 
tromagnetic coupling constant for monopoles {am — fif^/^c = 137/4n^ for a 
Dirac monopole). Nevertheless, the scattering cross sections for monopoles inci- 
dent on charged particles have been computed using non-perturbative techniques 
(see ref. [5] for a review of monopole scattering calculations). This has allowed 
the calculation of the stopping power of monopoles in material [246-249]. In addi- 
tion the binding energy of monopoles to nuclear magnetic dipoles has been calcu- 
lated [250, 251]. Since these calculations concern the interactions of monopoles in 
a detector, they are discussed in 5.1.2 and 6.4.2, respectively. 




Fig. 9. A Feynman diagram of a multi-photon process mediated by an internal monopole 
line in pp collisions. 
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Even though a perturbative treatment of monopoles is impossible, experiments nev- 
ertheless must assume certain production processes in order to estimate the detec- 
tor acceptance. The production mechanisms which have been assumed by experi- 
ments are typically those which would correspond to leading order pair production 
in a perturbative picture. A Drell-Yan-like mechanism to produce a monopole- 
antimonopole pair is a commonly considered scenario in hadron-hadron collisions 
[252]. Searches at e+e~ experiments have hypothesised monopole production through 
the annihilation reaction e+ + e~ — > Z'^/l — ^ [253]. In lepton-hadron 
collisions, pair production via photon-photon fusion has been envisaged [254]. 
To obtain calculations of cross-sections and of the kinematic distributions of pro- 
duced monopoles, the formalism for the electroweak production of is used, 
with appropriate phase space modifications to account for the larger monopole 
mass. However, it is important to emphasise that little is known about processes in 
which monopoles could be directly produced, and that the above reactions are 'best 
guesses'. Thus, some experiments also use random phase-space models to calculate 
acceptance. Examples of the methods used by experiments to model monopole pair 
production are discussed in Section 7.2. 

It has also been proposed in Refs. [255] and [256], respectively, that virtual monopoles 
mediating e+e~ and hadron-hadron collisions would be manifest via photonic ra- 
diation, as illustrated in Fig. 9 for pp collisions. Calculations of cross sections of 
such process assume that perturbation theory can be used for monopoles. Lim- 
its [257, 258] obtained using these calculations, which are discussed later in Sec- 
tion 7.2.5 have been criticised [259]. 



5 Interactions of SMPs 

This section describes the expected interactions of SMPs in particle detectors. The 
dominant types of interactions are electromagnetic and strong interactions. An out- 
line is given of the theory of electromagnetic energy loss for electrically and mag- 
netically charged particles. Compared to the electromagnetic case, nuclear interac- 
tions of SMPs containing heavy coloured objects are poorly understood. A review 
is given here of various models which have been proposed to describe these nuclear 
interactions and their uncertainties are also pointed out. 

5.1 lonisation energy loss 

The most commonly used observable in SMP searches is the measurement of the 
continuous ionisation energy loss Both electrically and magnetically charged 
SMPs lose energy principally through ionisation energy loss as they propagate 
through matter and for both types of particle the theory of electromagnetic energy 
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loss is well established. 



5.1.1 lonisation energy loss by electrically charged particles in material 

As an electrically charged particle moves through a material it loses energy either 
by interactions with atomic electrons or by collisions with atomic nuclei in the 
material. The first of these results in the liberation of electrons from the atoms in the 
material (i.e. ionisation) while the second results in the displacement of atoms from 
the lattice. The energy loss due to the second process is called the Non-Ionising 
Energy Loss (NIEL). The differential energy loss (stopping power), (dE/dx), due 
to the ionisation energy is much larger than the NIEL [260] in practical particle 
detectors. 

For fast particles of charge Zi in a medium of atomic number Z2 the mean ionisa- 
tion energy loss is given by the Bethe-Bloch formula [1]. 



dE 47re^Zf /I f2me&j3'YTn 



0-/5'- I (18) 



dx rriec'^P^ V ^ 



where e and me are the charge and mass of the electron, n is the number of electrons 
per unit volume in the material, (3 is the relativistic velocity of the incident particle, 
7=1/ Vl — and le is the mean ionisation potential of the material. The latter 
can be parameterised by [261]; 

Ie{Z2) = (I2Z2 + 7) for Z2 < 13 or (9.76^2 + 58.8^2^°-^^) eV for Z2 > 13. 
The quantity T„iax is the maximum kinetic energy which can be imparted to a free 
electron in a single collision and is given, for a particle of mass M, by 

l + 27me/M + (me/M)2" ^ 

The term S represents the density effect which limits the relativistic rise at high 
energy and it has been calculated by Stemheimer et al. [262]. This term is only 
relevant for particles with 3 and for massive particles Tmax ~ 2meC^/3^7^. 

Slight differences occur for positive and negative particles moving with low veloc- 
ity [263]. 

For low energies when the velocity of the incident particle is comparable or less 

than the velocity of the electrons in the atom, the so-called Lindhard region, this 
formula is no longer valid. The energy loss is then proportional to the particle ve- 
locity (3 [264]: 

dE 2 Z,Z2(3 

— = TV^eSTre a,,^— — , (20) 
dx Z Po 

where N is the number of atoms per unit volume, ~ zl^^, Qg is the Bohr Radius 

222 2 

of the hydrogen atom and Za — Zl + Z^ . This formula holds for 13 < Zl 13 o where 
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Po = 2^7^ (~ .0073) is the electron velocity in the classical lowest Bohr orbit of 
the hydrogen atom. 

The intermediate region, in which neither the Bethe-Bloch formula (Eq. 18) nor the 
Lindhard formula (Eq. 20) are valid, is defined by the velocity range Pi < P < P2, 
where pi = max[aZyV(l + «^!^^), (2^2'^ + l)/400], p2 = aZi/{l + aZi) and 
a is the fine structure constant. This region is described by Anderson and Barkas 
[265]. In this region a polynomial can also be used to join up the two regions [261] 
of the form; 

HE 

— = AP'' + Bp^ + Cp + D (21) 

ax 

where A, B, C, D are derived from the four simultaneous equations obtained by 

equating ^ and ^y^' at Pi and P2 for the Eqs. 18, 20 and 21. This gives: 

A - -J— (mi. k - '^^A 
^ ~ A/32 yi3^ -r rv J 

B = +P2)^p- iPi + 2/^2) f - Ay'k) 

C = ^{piiPi + 2P2)k + _ 

D = AFfe('^1^2(3/?2 - Pi) + Plyi{P2 - m) - PlP2{Plk + P2fX 

where A/? = P2 — Pi, y 1,1/2 are the ^ values computed from Eqs. 18 and 20 at ve- 
locities /3i and /32, respectively, Ay = yi-y2, k ^ [2^^^77,(7! -In ( ^7(i-^|)'' ))] 
and 72 = 1/^1 -/?|. 

The variation of ^ with P'y is illustrated in Fig. 10 by a calculation made for /i+ in 
copper [1]. The rise with P^ in the Lindhard region from Eq. 20 turns into a fall as 
P^ increases, as expected from Eq. 18. The accuracy of the rather arbitrary polyno- 
mial procedure between the two regions can be assessed from extrapolating the rise 
at low P'j and the fall at higher p-j into the worst possible case of a discontinuous 
join. The variation of ^ with material is shown in Fig. 1 1 [1]. 



5.1.2 lonisation energy loss by magnetically charged particles in material 

Since a Dirac monopole has a magnetic charge equivalent to an electric charge 
of 137e/2, such a particle would be expected to suffer ionisation energy loss at 
a rate several thousand times greater than that of a particle with electric charge e 
[246,247]. A Dirac monopole would thus be expected to leave a striking ionisation 
signature. lonisation energy loss could also lead to a monopole becoming stopped 
in detector material before reaching a tracking chamber. 

The energy loss in the form of ionisation of magnetic monopoles passing through 
material due to interactions with atomic electrons has been shown to follow a form 
similar to the Bethe-Bloch equation (Eq. 18) but without the multiplicative factor 
1/ P"^ [246,247]. The velocity-dependent Lorenz force is responsible for the inter- 
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Fig. 10. Stopping power (^) for positively charged muons in copper as a function of 
/?7 = p/M taken from the Particle Data Group [1]. The solid curve indicates the total 
stopping power. The radiative effects apparent at very high energies are much less relevant 
for particles heavier than muons. The different regions indicated by the vertical bands are 
described in Ref. [1] as are the small difference between positive and negative charges at 
low values of (the Barkas effect [266]) shown by the short dotted lines labelled /i^. 

action between a moving magnetic monopole and an atomic electron in the mate- 
rial. In contrast, for a moving electrically charged particle the velocity-independent 
Coulomb force is responsible for the interactions. The velocity dependence of the 
Lorenz force causes the cancellation of the 1 factor in the formula for ^ for 
magnetic monopoles (compare Eqs. 18 and 22). The detailed formula for the stop- 
ping power of a magnetic monopole of strength g is given in [246] as 

= -Inf ) + ^'^'^ -B{\g\) (22) 

and the modifications at very low velocity in [247]. Here is the mean ionisa- 
tion potential for magnetic monopoles which is close in value to Jg (Eq. 18). The 
relationship between and can be expressed as Im = /eexp— D/2. Stern- 
heimer [267] has shown, for several solids, that Dill) = 0.34, D{C) = 0.22, 
D{A\) = 0.056, D{Fe) = 0.14, D{Cu) = 0.13, and D{W) = 0.07. The correc- 
tion terms K{\g\) = 0.406, 0.346 and B{\g\) = 0.248, 0.672 for = 1, 2 Dirac 
Monopole strengths, respectively [246]. Fig. 12 (left) shows the stopping power 
for a unit Dirac magnetic monopole in aluminium as a function of the velocity of 
the monopole. Inspection of Eqs. 18 and 22 shows that the ratio of the stopping 
power for a unit Dirac monopole and a unit electric charge moving with velocity (3 
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Fig. 11. The stopping power (^) for particles in different materials (omitting radiative 
effects) [1]. 

is ~ 4700/3^. It can be seen from Fig. 12 that, as a monopole slows down, the ion- 
isation becomes less dense, in contrast to the case of electrically charged particles 
for which the opposite is true. This adds to the striking nature of a track left by a 
monopole. 

The large differential ionisation energy loss of monopoles makes it relevant to dis- 
cuss the range. The range R, of monopoles in aluminium is computed [254] by 
integrating the stopping power shown in fig 12. So that 



Fig. 12 (right) shows the computed range (normalised to mass), for a Dirac monopole 
versus P/M = (5^ where P and M are the momentum and mass of the monopole, 
respectively. 

The theory of energy loss described above is implemented as part of a GEANT 
package to describe monopole interactions in a detector [268]. 

5.2 Nuclear interactions ofSMPs 

As heavy (charged or neutral) hadrons propagate through a medium, they may un- 
dergo scattering from the nuclei of the material of the apparatus. For R hadrons 




(23) 
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Fig. 12. Left: the ^ for a Dirac monopole in aluminium as a function of the velocity of the 
monopole taken from [246] and adjusted for the electron density in aluminium. Right: the 
ratio of range to mass for a Dirac monopole in aluminium versus /?7, calculated from the 
stopping power, ^. 

the interaction cross sections are expected to be of the order of those for pion scat- 
tering (see below). As is shown in this section it is expected that the energy loss 
of an i?-hadron through such scattering can be sufficiently small as to allow it to 
penetrate through to an outer muon chamber and be reconstructed as a slow mov- 
ing exotic particle. However, in extreme cases such interactions can have a large 
impact on experimental searches. For example, these interactions can lead to R- 
hadrons being 'stopped' inside the detector material, because they come to the end 
of their ranges due to electromagnetic and hadronic energy losses. Also, R-hadrons 
could undergo charge exchange reactions in a hadron-absorbing material such as 
a calorimeter, in which case the initial charge at the production vertex is not nec- 
essarily the same after traversing the material. An understanding of such effects is 
crucial to quantify a discovery or assess a discovery potential. This section provides 
an overview of the theory of i?-hadron scattering processes and a description of the 
phenomenological approaches which have been used to describe them. Before dis- 
cussing these different approaches, we summarise which general observations can 
be made, independently of any specific model. 

One important feature of i?-hadron scattering common to all approaches is the pas- 
sive nature of the exotic heavy coloured object. The probability that the parton Ci 
of colour state i will interact perturbatively with the quarks in the target nucleon 
is small, since such interactions are suppressed by the squared inverse mass of the 
parton. As a consequence, the heavy hadron can be seen as consisting of an essen- 
tially non-interacting heavy state Ci acting as spectator, accompanied by a coloured 
hadronic cloud of light constituents, responsible for the interaction. Hence the in- 
teraction cross section will be typical of that of a meson. In addition, the effective 
interaction energy of the heavy object is small. As an example, consider a Csqq 
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Fig. 13. i?-hadron-proton scattering processes, (a) Elastic scattering, (b) Inelastic scattering 
leading to baryon and charge exchange, (c) Inelastic scattering leading to charge exchange, 
(d) Resonance formation. 

state with a total energy E=450 GeV and a mass m of the Cg parton of 300 GeV, 
the Lorentz factor will be 7=1.5. Although the kinetic energy of the i?-hadron is 
150 GeV, the kinetic energy of the interacting qq system is only (7 — l)mqq 0.3 
GeV, (if the quark system consists of up and down quarks). For i?-hadrons produced 
at the Tevatron or LHC with masses above 100 GeV, the centre-of-mass energy of 
the system of quarks and a stationary nucleon can thus be at most around a few 
GeV. Thus, the energy scales relevant for heavy hadron scattering processes from 
nucleons are low and comparable with low-energy hadron-hadron scattering for 
which Regge theory is often applied. The heavy state Cj serves only as a reservoir 
of kinetic energy. 

Although i?-hadrons may scatter elastically or inelastically the energy absorbed in 
an elastic scattering process, such as that illustrated in Fig. 13 (a), is expected to be 
small [269], since the high-mass i?-hadron scatters on a lower mass target nucleus, 
and inelastic collisions are expected to be largely responsible for the energy loss 
of an /?-hadron. These inelastic collisions may cause the conversion of one species 
of /?-hadron to another in two ways: baryon exchange, which was overlooked until 
recently [235], and charge exchange, as shown in Fig. 13 (b) and (c), respectively. 
In the first process, an exothermic inelastic i?-meson-nucleon interaction results in 
the release of a pion. The reverse reaction is suppressed by phase space and because 
of the relative absence of pions in the nuclear environment. Thus, most i?-mesons 
will convert early in the scattering chain, in passing through hadron absorbing ma- 
terial, e.g. a calorimeter, to baryons and remain as baryons. This is important, since 
baryons have larger scattering cross sections. Baryon formation offers one oppor- 
tunity for a charge exchange process to take place. Charge exchange may arise 
in any meson-to-meson, meson-to-baryon, or baryon-to-baryon process. Although 
exact predictions of individual processes are difficult to make, the low energies in- 
volved in i?-hadron scattering imply that reggeon and not pomeron-exchange will 
dominate, and thus charge exchange reactions may well form a substantial contri- 
bution to all interactions. This may lead to striking topologies of segments of tracks 
of charged particles with opposite signs of charge on passage through hadron ab- 
sorbers or calorimeter material. It is also interesting to note that such a configuration 
can also arise if a neutral i?-meson, formed as an intermediate state during scatter- 
ing, oscillates into its own anti-particle and then subsequently interacts to become 
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a charged i?-hadron [244,245]. 



Several phenomenological approaches have been developed [29, 235, 270] to de- 
scribe _R-hadron nuclear scattering which are described later in Sections 5.2.1 and 
5.2.2. Although these differ in the phenomenology used, they are largely based on 
the generic picture of i?-hadron scattering described above, much of which was 
first introduced in [269]. Low-energy hadron-hadron data are typically used to es- 
timate scattering cross sections and several of these models are based on modified 
Regge paramaterisations of the data. Uncertainties in the models arise from several 
sources. Most models assume _R-meson scattering processes to be dominant despite 
it being likely that a meson will convert into a baryon and then stay baryonic as it 
propagates through matter. A further theoretical uncertainty arises from resonance 
production. The formation of resonant i?-hadron states, as shown in Fig. 13 (d), is 
expected to take place during hadronic interactions. To date, no explicit modelling 
of i?-hadron resonances has yet been attempted. However, it has been argued that 
the the minimum centre-of-mass energy required for a scattering process is above 
the expected mass of the main resonances [235]. Nevertheless, since little is known 
about such resonances it is important to note that they may play an important role 
in accounting for i?-hadron energy loss and conversions. 

Below, we will discuss for each of the phenomenological approaches the interaction 
cross section or interaction length, scattering processes included, and energy losses. 
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Fig. 14. Predictions from four phenomenological models of expected hadronic energy loss 
per interaction as a function of the Lorentz factor 7. Also shown is the ionisation energy 
loss corresponding to the passage of an i2-hadron with charge ite through 18 cm of iron. 
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5.2.1 Models based on Regge theory 



Consider a gluino i?-hadron scattering proces RN — > R'X, where R is the initial 
R-hadron before scattering, N isa single nucleon (proton or neutron), R' is the R- 
hadron after scattering, and X is the remaining system produced in the scattering, 
the latter is expected to be a nucleon and possibly a number of other light hadrons. 
Baur, Cheung and Gunion (BCG) [270] propose three different ansatzes for the 
functional form of the cross section for gluino i?-hadron scattering processes: 

(1) ^^i^ oc 1 for \t\ < 1 GeV^ and for \t\ > 1 GeV^ 

(2) ^|^|^2 given by a triple-pomeron form used to describe single inclusive par- 
ticle production in pion-nucleon scattering 

^ cx ^/?2(|i|)(^)2("^(l*l)-i)(m^)"^(o)-i where ap{\t\) = 1 - 0.3|t| 



da 



d\t\dm^ 

and/?(|t|)- ' 



(1+1*1/0.5)2- 

(3) ^^i^ oc 1 for \t\ < 4 GeV2 and for \t\ > 4, 

where t is the four-momentum transferred from Rto R', and nix is the mass of the 
remaining system of final state particles produced. The cross section axiRN) for 
gluino i?-hadron scattering on a nucleon is derived from the cross section axiTrN) 
for pion nucleon scattering. The normalisation of the cross section is determined by 
constraints on the collision length Xt{R) of an i?-hadron. The evaluation of Xt{R) 
is made by correcting the pion collision length AT(7r). 



i 



(ZIL)XjlL (24) 



The colour factors Cp — 4/3 and Ca — ^ arise due to the low-mass colour-octet 
constituent of the i?-hadron. The term (r^) is the squared transverse size of the 
particle. In the case of a pion and an /?-hadron is given by (r|,) oc l/m^g, and 
{rl) oc Here, rriq and rrig are the constituent masses of light quarks and gluon 

constituent masses, respectively, which were assumed equal in this approach. This 
leads to a collision in length in iron of i?-mesons of around 19 cm. However, the 
uncertainties in the constituent masses of the partons would affect the collision 
length, and a value of 38 cm has also been considered within the BCG approach. 

As mentioned above, this model does not take into account i?-baryons or conver- 
sion from mesons into baryons. Charge exchange reactions are possible since the 
i?-hadron is considered to be stripped of its system of quarks following a nuclear 
scattering. A refragmentation process, governed by a probability to fragment into a 
charged or neutral state then allows the formation of a new type of i?-hadron. 

The energy loss of hadron it! scattering on nucleon N in the process RN — > R'X 
is given by 

AE = <-^N + \t\_ (25) 
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Here AE is then evaluated according to the parameterisation chosen for the differ- 
ential cross section ^pj^;^- The distribution of mean energy loss per collision as a 
function of 7 for the three different scenarios of this approach is shown in Fig. 14. 
All three approaches show a rising energy loss with 7. However, there are large 
differences between the various scenarios. For comparison, the ionisation energy 
loss corresponding to the passage of an it!-hadron with charge ±e through 18 cm of 
iron (1 interaction length for R-mesons in the approach by Kraan discussed below 
in Section 5.2.2) is also shown. 

This approach has been extended by Mafi and Raby [29] (MR), where again a sin- 
gle particle inclusive scattering process RN R'X involving only i?-mesons was 
considered. However, MR consider two Regge trajectories: an isosinglet pomeron, 
and an isovector reggeon p trajectory. Using these, triple pomeron and reggeon- 
reggeon-pomeron cross-section forms were extracted. The presence of reggeon ex- 
change incorporates charge exchange processes naturally without relying on the 
BCG assumption that an 7?-hadron is stripped of its soft partonic system in an in- 
teraction and then forced to fragment. Two values of the collision length A^- = 19 
cm and A^- = 38 cm are used in their analysis to estimate i?-hadron stopping in 
different scenarios and the scattering cross section functional form was rescaled to 
achieve these values. The relative proportion of reggeon and pomeron were also 
allowed to vary between the extreme cases of 100% reggeon and 100% pomeron 
exchange. As mentioned above, charge exchange processes are naturally included 
in this approach. Baryon exchange processes are however omitted. The typical en- 
ergy loss per collision is of the order of several GeV and is comparable with the 
BCG approaches. 



5.2.2 Model based on geometrical cross sections 

Since the behaviour of the hadron-hadron scattering cross sections at values of 
the centre-of-mass energy below several GeV is specific to the type of hadrons 

interacting, it is not necessarily a reliable general guide to the scattering of R- 
hadrons. Therefore, in a complementary approach to the Regge-based models of 
MR and BCG the constant geometrical cross section is used by Kraan [235] over 
the full scattering energy regime. 

The total nucleon interaction cross section is approximated by the asymptotic val- 
ues for the cross sections for normal hadrons scattering off nucleons. The model 
assumes that only u and d quarks are present in i?-hadrons and that each quark 
which can interact represents a contribution of 12 mb to the total scattering cross 
section. Thus, the scattering cross sections of a gluino i?-meson and i?-baryon are 
24 mb and 36 mb, respectively. A gluino-gluon state can be assumed to have the 
same cross section as a gluino i?-meson, since the geometrical cross section is ap- 
proximated by the high-energy hadron cross section, where gluon exchange would 
dominate. The gluon-gluon coupling is a factor 9/4 larger than the quark-gluon 
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coupling, but a meson has two quarks, resulting in a cross section of a gluino-gluon 
state which is (9/4)/(l+l):^l times the cross section for a gluino R-meson. Trans- 
lating this into interaction lengths, for i?-baryons the average nuclear interaction 
length (i.e. amount of material where on average one interaction takes place) is 12 
cm in iron, 3 1 cm in carbon, and 660 cm in hydrogen. For i?-mesons these numbers 
are 3/2 larger. 

This model includes predictions for all possible 2^2 and 2^3 processes. The 
relative rates of 2 — > 2 and 2^3 processes are, in the asymptotic region, set to 
15% and 85%, respectively, as suggested by hadron-hadron scattering data. A phase 
space factor is used to determine the proportion of 2 — > 2 and 2 — > 3 processes 
as the scattering centre-of-mass energy is reduced. Within the sets of 2 2 or 
2^3 interactions each allowed process is assigned the same matrix element and 
the different rates of processes is determined by phase space. Included in these 
processes is charge exchange, and, for the first time, baryon exchange. Fig 15 shows 
the fraction of _R-mesons which convert into i?-baryons when travelling in iron. As 
was already mentioned in Section 5.2, conversion from baryons back to mesons is 
highly suppressed. 

Energy losses are determined by implementation of the model into the GEANT 
framework (GEANT-3 [235] and GEANT-4 [271]), which allows a more sophisti- 
cated treatment of energy loss in a nuclear reaction than was possible in the ana- 
lytical approaches of BCG and MR. Issues related to nucleus scattering like Fermi 
motion, binding energy of the nucleons inside the nucleus, evaporation energy and 
instability of a nucleus are thus included. Figure 14 displays the nuclear energy loss 
per interaction for an i?-meson (the curve labelled Kraan). The energy losses from 
this model are found to lie between the different BCG approaches. lonisation losses 
are also shown and they dominate energy losses only at low beta values (/? < 0.75). 
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Fig. 15. Fraction of i?-mesons converted into i?-baryons as predicted in the model of 
Kraan [235]. 
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6 Search techniques 



Following the previous section on the interactions of SMPs in materials, this section 
describes how these interactions can be used to search for SMPs in colliders. Search 
techniques include the use of ionisation energy loss, Cherenkov radiation and time 
of flight. Finally, we summarise search techniques which are only applicable to 
magnetically charged SMPs. 



6.1 Methods based on ionisation energy losses 



The most commonly used observable in SMP searches is the measurement of the 
continuous ionisation energy loss 



6.1.1 Measurements of ionisation energy losses in tracking systems 

Ionisation energy loss can be used to identify heavy charged particles with track- 
ing chambers. The measurement of ^ is part of the routine program of calibrating 
charged particle tracking chambers. However, a number of systematic studies must 
be performed in order to optimise the detector calibration. Examples of the appli- 
cation of ionisation loss measurements in different collider environments can be 
found in Section 7.1, and the accompanying table 5. 

It can be seen from Eq. 18 that the value of ^ is dependent on the velocity fac- 
tors (3, 7. The particle's momentum {p = M(3'^) is measured independently, usu- 
ally from the curvature of the track in the magnetic field. Hence, comparing the 
measured value of ^ with an independent measurement of track momentum, the 
particle mass can be determined. Fig. 16 shows an example of what mass separa- 
tion was achieved in the OPAL experiment when simultaneously measuring ^ and 
momentum. The regions in which SMPs possessing various values of charge and 
mass could be manifest are shaded. 

Several issues must be taken into account in the identification of charged particles 
by 4^ measurements. First of all, considerable fluctuations occur in single measure- 
ments of These are not Gaussian but are asymmetric with a high-energy tail due 
to the emission of energetic secondary electrons {5 rays), giving rise to the Landau 
distribution [272]. Examples of methods adopted to avoid complications due to 
these high-energy tails include disregarding abnormally large single measurements 
of ^ and taking the mean of the remaining measurements [273] or suppressing tail 
distributions [274]. Another method is to perform a maximum likelihood fit of a 
Landau distribution to the sample of measurements on a track [275] to determine 
the mean value of 4^ for the track. The particle masses determined from the 4^ 

ax ^ ax 

at a fixed momentum had smaller high-mass tails using this method than with the 
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others [275]. Hence this is likely to be the best technique to employ in searches for 
unknown particles, since there will be less background by smearing from the lower 
mass region. The large numbers of measurements needed for reliable ^ determi- 
nation means that silicon vertex detectors with small numbers of detector layers are 
likely to be of limited usefulness in searches for unknown SMPs. 

Further problems in measuring ^ for unknown massive particles could arise due 
to saturation of the electronics used in the track detectors. For such searches the 
electronics must have a wide enough dynamic range otherwise large values of 
^ become unmeasurable due to saturation. This is particularly true for magnetic 
monopoles, for which ^ is significantly larger than that for electrically charged 
particles, but may already play a role for slow electrically charged ones. Although 
saturated hits have a tracking resolution much worse than that of typical unsatu- 
rated data, the presence of many saturated channels is itself a distinctive feature of 
the presence of a high-mass charged particle. A characteristic pattern of a charged 
particle helix can generally easily be recognised from the spatial distribution of 
saturated hits. The technique of searching for saturated hits has been applied in 
Ref. [276]. Saturation of electronics as a consequence of highly ionising particles 
may in particular have a non-negligible effect in the future LHC collider experi- 
ments where the resulting dead time as a result of electronics saturation may be of 
the order of the bunch crossing time. The effect of highly ionising particles on the 
CMS silicon strip tracker has been studied in Refs. [277, 278]. 

There are several kinds of backgrounds producing highly ionising signals. An im- 
portant background source to such studies is positively charged nuclei. These typi- 
cally arise as spallation products from secondary interactions of particles produced 
in the primary interaction. Such interactions take place with nuclei in the material of 
the apparatus. This accounts for stronger limits for negatively charged SMPs using 
this approach, since spallation products are positively charged. Another background 
arises from the finite resolution of a tracking chamber, causing, for example, two 
overlapping tracks to be measured as one; this results in a highly ionizing signal. 
These backgrounds are mainly caused by photon conversions. Finally, interactions 
of the colliding beams with residual gas atoms in the beam pipe can also produce 
highly ionising spallation products. However, most background events can usually 
be effectively removed by the application of a lower cut on the momentum [279]. 



6.1.2 Methods using ionisation damage in plastic foils 

Next to ^ measurements to identify heavy particles in tracking chambers, heavily 
ionising single particles (magnetic or electrically charged) can be detected by the 

damage they cause to certain plastic foils eg. foils made from plastics such as lexan, 
CR39 and Makrofol [127]. Such damage is caused by both the NIEL and the ioni- 
sation energy loss. The effect of the damage is made visible by chemically soaking 
the foil in a concentrated alkali solution which etches the region around the damage 
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Fig. 16. Measurements of ^ vs ^ as measured by the OPAL experiment [280], where Q is 
the charge of the particle. The smooth curves show the expected values of the mean values 
of ^ for the known particles, while the dotted points are the measured values. The regions 
in which unknown massive particles are sought are shaded. Also shown are lines which 
illustrate the expected ^ values of SMPs with specific masses and charges. 

centre into small holes (pits) which are visible under a microscope [281]. The value 
of ^ of a heavily ionising track can be determined from the etch rate i.e. the ratio 
of the rate of increase in depth of the pit to the decrease in overall thickness of the 
foil due to the chemical action. Fig. 17 shows the pits produced by sulphur ions 
in a heavily etched CR39 foil (type ENS) as reported in [282]. The etching pro- 
cess only makes the tracks visible for ionisation levels above a threshold value so 
that the lightly ionising tracks remain invisible and only the heavily ionising tracks 
are seen. For heavy ions of velocity (3 and atomic number Z, the experiment of 
Ref. [282] found that the threshold corresponded to ions with Z / (3 S for CR39 
(type ENS). In collider experiments [25S, 282-286], the layers of foils surround 
the beam pipe at an interaction point and are left to be exposed to the products of 
the beam interactions. After the exposure the foils are soaked in the appropriate 
chemical to etch out the damage centres. Usually, the inner layer foil is scanned 
under a microscope to search for etched pits. The outer layers are then scanned for 
pits which align with those found in the inner layer to search for heavily ionising 
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continuous tracks. 




Fig. 17. Microphotograph of a heavily etched CR39 sheet (type EN3) exposed to sulphur 
ions of 200 GeV/nucleon [282]. The sheet thickness was reduced from 1.4 mm to 0.2 mm 
by the heavy etching. Note the holes produced by the sulphur nuclei. 



6.2 Ring Imaging Cherenkov detection 

Another way to exploit the electromagnetic interactions of SMPs to identify them 
is via the Cherenkov effect [1,287]. A particle, travelling with velocity [3 which is 
faster than the speed of light in a medium of refractive index n, emits light waves 
by the Cherenkov effect in a cone at an angle of cos 9 — l/n/3 relative to the track 
direction. The number of Cherenkov photons per unit path length is proportional 
to sin^ 9. Particles with velocity less than the speed of light in the medium do not 
emit Cherenkov light. Unlike the other LEP experiments, the DELPHI detector 
contained Ring Imaging Cherenkov (RICH) systems [288-291] comprising two 
Cherenkov radiators with different refractive indices: one in the liquid and one in 
the gaseous phase. To select low-velocity massive particles it was required that 
neither detector should give a detected Cherenkov light signal for high-momentum 
tracks. The DELPHI experiment was able to use this technique for a number of 
SMP searches [292-294], the results of which are described in Section 7.1. 

6.3 Methods based on time of flight 

Another technique which can be used for the detection of SMPs is the method of 
time of flight. Although less widely used than the ^ method, this method has been 

applied in several collider environments, examples of which can be found in Sec- 
tion 7. Massive particles are produced with a smaller velocity than light particles, 
and thus would have a larger time of flight. For electrically charged particles the 
mass of the particle can be determined by correlating the flight time with the track 
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momentum, p, measured from another source e.g. from the track curvature in a 
magnetic field. Timing can be obtained from the detection of the tracks in scintilla- 
tion counters. Since M = p/ j3'^, particles are identified by plotting a graph of 
against l/p^ which is linear with slope for particles of mass M. The accuracy 
of determination of the mass due to the momentum and time of flight uncertainties, 
which are assumed to be uncorrelated, is 



where the uncertainty in the velocity Ap/P = At/t with t the time of flight. For 
a typical timing resolution of At ~ 1 ns and performing time of flight over a 
distance of ~ 3 m [295], the mass resolution AM/M varies between 2% to 23% 
for particles in the range of 0.2 < (3 < 0.8, neglecting the momentum resolution. 
Hence the technique is comparable to the ^ technique in the precision of the mass 
measurement [275]. 

Backgrounds in time of flight measurements can arise from instrumental effects 
such as mismeasured times or random hits in the scintillation counters. Cosmic 
rays may form a substantial background, as may i^-decays in flight. In addition, at 
colliders with high-frequency bunch crossings and a large number of interactions 
per bunch, such as the LHC, the particles from one bunch may become confused 
with particles from another bunch. In this situation, it may not be possible to make 
a unique determination of the time of flight. 

A slightly different technique which has not been used so far but could be used 
to detect heavy charged particles is to measure the track velocity with a tracking 
system based on wires and gas. The detection time of the signal of a traversing 
charged SMP is in simplified form given by 



where to is bunch crossing time, troF is the time of flight from the beampipe to the 
wire, tdrift is the drift time of the signal inside the gas to the wire, and teiectronics is 
the time for the transmission of the signal. For a slow particle, tTOF is large. Since 
default track reconstruction programs are based on expectations for txoF for light 
particles, a misalignment occurs for slow moving particles. If the reconstruction of 
the track is successful, this pattern is very distinctive for heavy slow particles. It 
could in addition be used to support a heavily ionizing signal. This method has not 
yet been applied in data, but has been studied in Ref. [296] with the ATLAS muon 
reconstruction software [297]. 




(26) 



^ — ^0 + trOF + tdrift + ^electronics 



(27) 
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6.4 Specific techniques for magnetic monopoles 

In addition to the techniques described above, a number of approaches have been 
used purely to search for particles with magnetic charge. Examples of applications 
of the techniques discussed below are given in Section 7.2. 

6.4.1 Parabolic tracks 

Monopoles experience a force giB in a magnetic field B which causes them to 
accelerate. With the field aligned along the z axis they assume a parabolic trajectory 
with 

^ ^|B|r^ r ^ ^g20.54|B|r^ r 

2ePT/5Tl09 tan^o '^Prh tan^o 

where is the z coordinate of the vertex and z{r) is the coordinate of a point on 
the trajectory at distance r from the proton beam with lengths in metres and |B| 
in Tesla. The transverse momentum and transverse velocity of the monopole are 
Pt GeV/c and (3t, respectively, and (^d is the pole strength in units of the Dirac 
Monopole. The initial angle of the monopole to the magnetic field direction is 9q 
and e is the unit of electric charge. In this equation g is the magnetic pole strength 
which is negative (positive) for South (North) poles which decelerate (accelerate) 
in the -\-z direction in the magnetic field. 

A simulation of the passage of a monopole, anti-monopole (mm) pair at the TASSO 
experiment [298] is shown in Fig. 18. Fig. 18 (top) shows the r — plane (transverse 
to the magnetic field direction) of the paths of the two particles which are produced 
at the primary interaction point (IP). However, as seen in the lower plots, the paths 
of the particles are clearly disturbed by the magnetic field in the s — z view, where 
2; is a distance parallel to the beam axis and s is the total distance travelled by the 
monopole along its path. 

6. 4. 2 Searching for stopped magnetic monopoles 

A further method to search for monopoles utilises the fact that the ^ of magnetic 
monopoles is so large that they tend to stop in the beam pipe or apparatus surround- 
ing the interaction region. When a monopole stops, i.e. reaches a speed comparable 
to that of an electron in a Bohr orbit of the atom (/? ^ 0.01), it is expected to 
become bound to the nuclei of the atoms of the material. The binding energy is de- 
pendent on the magnetic dipole moment of the nucleus [299-301] which depends 
on the nuclear spin. Hence magnetic monopoles should remain bound in materi- 
als such as aluminium (nuclear spin 5/2). Since the binding energy of a magnetic 
monopole to a nucleus with zero magnetic moment is expected to be small, the 
monopoles may diffuse out of materials made of mainly even-A and even-Z nuclei 
with spin zero i.e. zero magnetic moment, e.g. carbon. The strong divergence of the 
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Fig. 18. Simulated monopole, anti-monopole production at the TASSO detector [298]. The 
top plot shows the paths of the particles (produced at the centre) in the r/(j) view. The 
lower left (right) plot shows the path of the monopole (anti-monopole) as it is accelerated 
(decelerated) in the magnetic field, in the s — 2; view, where s is the total distance travelled. 

monopole magnetic field, B = / ^Tir"^ , causes a persistent current to flow if the 
monopole is passed through one or more superconducting (sense) coils . In con- 
trast, divergenceless magnetic fields fi-om ubiquitous magnetic dipoles and higher 
moments cause the current to return to zero after complete passage through the coil. 
If the material surrounding the beam is cut into small pieces which are then passed 
through the superconducting coil a residual persistent current signifies the presence 
of a monopole in the sample. The current is measured using a Superconducting 
Quantum Interference Device (SQUID) connected to the sense coils. This method 
was invented to search for monopoles in lunar rocks [302-305]. The response of 
the superconducting coils can be calibrated by the traversal of the coil by one end 
of a long, thin solenoid, the magnetic field at the end of which approximates to 
that of a magnetic monopole. Since the induced monopole current is persistent, it 
is also possible to increase the sensitivity of the apparatus by repeated traversals 
of the same sample of material. This technique of gradually building up a signal is 
especially useful when searching for particles with values of the magnetic charge 
substantially less than the Dirac charge go- Using this technique experiments have 
searched for magnetic charges as low as [254,302]. The apparatus used for 
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Fig. 19. A schematic diagram showing the SQUID apparatus used by the HI experi- 
ment [254]. The conveyer belt travelled in small steps until the sample was passed com- 
pletely through the coil. At each step the current in the superconducting coil was read. 

the HI search [254] is typical of such searches and is illustrated in Fig. 19. 

In addition to using a SQUID, it is also possible to use strong external magnetic 
fields to look for stopped monopoles. In Refs. [306-308] materials in which monopoles 
could be stopped were placed in front of a pulsed solenoid. A magnetic field of 
around 80 kG is expected to be large enough to liberate a monopole and accelerate 
it towards the detection systems comprising scintillators and plastic track detectors. 
Searches have been sensitive to charges in the range 0.03-0.245(£). 



7 Searches at colliders 



This section gives an overview of the various direct searches which have been per- 
formed for SMPs at colliders, i.e. searches based on detecting SMPs by their pas- 
sage through the detector. The results of direct SMP searches broadly fall into three 
different classifications. 

• The first class concerns general searches for SMPs, which are made without 
assumptions on the properties of the SMP, other than the values of the electric 
or magnetic charge, spin and mass. Predicted cross sections are not available 
and model-independent upper limits on production cross sections are extracted 
within a specific mass acceptance region. Such searches are designed to be sen- 
sitive to SMPs with a wide range of values of mass and charge, and thus cover 
many of the SMPs predicted in Section 2, as well as unexpected particles. 

• The second class of searches assumes the existence of certain types of SMPs 
with specific quantum numbers. These searches are based on minimal theoret- 
ical scenarios in which the production mechanisms and mass -dependent cross 
sections are assumed to be known. The influence on the SMP production cross 
sections of any other exotic particles which may be predicted within these sce- 
narios, eg via loop diagrams, is small. This allows lower limits on SMP masses 
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to be derived which depend dominantly on the properties of the SMP and not on 
other model parameters. 
• The third class concerns searches made within theoretical scenarios for which 
the search results are quoted in terms of the complicated model parameter space. 

Searches for electrically charged particles have been made in all three classes whereas, 
owing to the difficulties in calculating short-distance monopole interactions, monopole 
searches are typically made within the first class, with comparatively few searches 
also quoting a mass limit. It should be noted that a given study can present sets of 
results which fall into different classes. 

This section is organised as follows. In Section 7.1 direct searches for electri- 
cally charged particles are reviewed. Searches at e+e~, hadron-hadron, and lepton- 
hadron facilities are described. Searches in each collision environment are de- 
scribed in the class order given above. We then make a short summary of all of 
the direct collider searches, pointing out the most stringent limits which have been 
obtained. Finally, some indirect searches, i.e. searches for physics signals in which 
SMPs play an indirect role (for example from their presence in loop diagrams) but 
may not manifest themselves as final-state particles, are briefly touched upon at the 
end of Section 7.1. The same structure is adopted for the description of monopole 
searches in Section 7.2, although in this case the vast majority of search results fall 
within the first classification. 

The intention in this section is not to provide an exhaustive compendium of all 
results in this field, such as that which can be found in Ref. [1]. However, tables of 
selected results showing the most stringent limits extracted for electrically charged 
particle searches which belong to class 1 and to classes 2 and 3 are summarised in 
Tabs. 4, 5, respectively. Tab. 6 summarises magnetic monopole searches. 

7. 1 Searches for electrically charged SMPs 

Searches for SMPs have been performed at e+e^, lepton-hadron, and hadron-hadron 
colliders. The majority of recent searches has been performed at LEP and compar- 
atively few studies have been made at the Tevatron and HERA. 

7.1.1 Searches at e+e" experiments 

The early studies at low-energy e+e" facilities in the 1980s, such as PETRA (a/s — 
27 - 35 GeV), PEP (^i = 29 GeV) and Tristan (^i = 50 - 61 GeV), com- 
prised generic searches for SMPs with unexpected mass and charge, and, in par- 
ticular, searches for particles possessing fractional charges, as inspired by the pos- 
sibility of the existence of free quarks [1]. Although the notion of the existence 
of free quarks had fallen out of favour in the 1990's, the experiments at LEP 
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(v^ = 91.2, 130 — 209 GeV) continued to search for fractionally charged objects, 
and, in addition, made a number of SMP searches within exotic scenarios, such as 
SUSY. Owing to the high collision energies and large luminosity samples at LEP-1 
(LEP-2), typically ~ 100 (~ 700) pb~\ the most stringent results were usually 
extracted at LEP. Unless stated to the contrary, the LEP results quoted here were 
based on luminosity samples of approximately these magnitudes. 

The e+e~ searches for generic SMPs used simple topological cuts to look for the 
exclusive or inclusive production of SMPs. The principal experimental observable 
used was the ionisation energy loss measured in a tracking system. This was the 
approach adopted by the low energy experiments [309-316] and by ALEPH [276, 
279], OPAL [280,317-319] and L3 [320-322]. For the DELPHI studies, a RICH 
detector was used, in combination with ionisation measurements from a TPC [292, 
293,323]. 

The most stringent e~^e~ limit, in terms of mass reach, for fractionally charged 
particles, was obtained recently at LEP-2 by OPAL [280]. For this study, the col- 
lision centre-of-mass energy spanned 130-209 GeV. An exclusive pair-production 
mechanism e+e" — > XX was assumed, and model-independent upper limits on 
the production cross section for weakly interacting scalar and spin | particles with 
charge ±e of between 0.005 and 0.03 pb were extracted at 95% confidence level 
(CL), for a sensitive mass region between 45 and 103 GeV. This work also made a 
search for particles possessing fractional charges ±|e, ±|e, ±|e, with the resultant 
cross-section limits lying between 0.005 and 0.02 pb at 95% CL. Using a substan- 
tially smaller luminosity sample (~ 90 pb~^), and for centre-of-mass energies of 
130-183 GeV, DELPHI was able to place limits on the production cross-section of 
SMPs with charges ±|e [293], for masses in the range 2-91 GeV. In this work a 
free squark model was assumed, and cross-section limits of between 0.04 and 0.6 
pb were obtained at 95% CL. Using an exclusive slepton pair-production model, 
DELPHI also extracted upper cross-section limits in the range 0.05-0.3 pb for SMPs 
with charge ±e for masses up to 93 GeV [292,293]. Less stringent cross section 
limits for SMPs with charge ±e were also extracted by ALEPH [276]. The L3 
experiment [320-322] searched for signatures of charged heavy leptons at LEP2. 
However, the results of these studies did not include any upper cross-section limits 
which are relevant for generic SMPs. 

Searches at LEP-1 for fractionally charged objects made by ALEPH [279,318] and 
lower energy e+e^ colliders typically expressed their results in terms of limits on 
Rx, the ratio of the cross section for the single or pair production of SMPs to that of 
exclusive production. Fig. 20 shows limits on Rx (to 90% CL), which were 
obtained from e+e" experiments for fractionally charged SMPs with masses below 
45 GeV [279,309-316] for the putative charges ±|e, ±|e, and ±|e. Also shown 
is the ratio of the upper limit on the single inclusive production cross section of 
an SMP with charge ±e to the exclusive dimuon cross section. The LEP-1 results 
from ALEPH [279] and OPAL [318] were based on luminosity samples of 8 and 74 
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Fig. 20. The upper Umits at 90% CL on the ratio Rx as extracted at LEP-1 (ALEPH [279] 
and OPAL [318]) and lower-energy experiments [309-316]. 

pb~^, respectively and extend down to 5-10 GeV in mass. The lower energy e^e' 
limits exclude SMPs with masses as low as ~ 1 GeV. To obtain the results shown 
in Fig. 20 the experiments assumed that the SMPs follow a momentum distribution 
suggested for a heavy particle: =constant [222] (see Section 4.3). The ex- 

tracted limits are very sensitive to the form of the momentum distribution. ALEPH 
demonstrated that the limits can change by more than a factor of five if the mo- 
mentum dependence is derived from a fit to the Feynman-x spectra of inclusively 
produced hadrons [279]. 



In addition to the results described above, OPAL also extracted mass-dependent 
cross-section limits on the single inclusive production of SMPs with charges ±e,± |e, 
±|e, and ±2e. However, the experiment was not sensitive to particles with charge 
±|e. Furthermore, ALEPH and OPAL and some of the lower energy experiments [309, 
3 10, 3 12, 3 14] also presented limits on the exclusive pair-production of SMPs with 
specific charges. The obtained limits on SMPs with charges ±|e, ±e, ±|e, 2e were 
of similar values to the limits they quoted for the inclusive cases. 
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The possibility of the production of SMPs with values of charge substantially in 
excess of e, such as Q-balls, has been investigated by experiments using plastic 
track detectors. While these were primarily monopole searches, the experiments 
were also sensitive to non-magnetically charged particles with unexpected values 
of electric charge. These works are often overlooked in studies which survey SMP 
searches although they are extremely important in extending the range in charge 
to which collider experiments are sensitive. Several such experiments took place 
at LEP-1 [253,286]. In Ref. [286] the MODAL passive detector was deployed at 
the vacant 15 intersection point at LEP and was exposed to e+e" collisions corre- 
sponding to an integrated luminosity of 48 nb~^. The experiment was sensitive to 
SMPs with electric charges of up to approximately |240e| at low masses, with de- 
creasing charge sensitivity as the mass increases. An upper limit on the production 
cross-section of highly ionising particles of 70 pb was extracted, assuming that the 
efficiency was equal to the maximum value of acceptance for a pair of highly ionis- 
ing particles. A second search [253], using passive detectors deployed at the OPAL 
intersection point, was exposed to a far higher luminosity (8.7 pb~^). However, 
this experiment quoted its sensitivity only for particles possessing various values of 
magnetic charge (see Section 7.2.1). 

In addition to the generic searches, model-dependent (class 2) searches for elec- 
tromagnetically charged and weakly interacting SMPs have also been made. Using 
the ^ technique heavy lepton limits have been extracted at LEP [280, 324] for a 
range of putative non-SUSY, heavy leptons, such as so-called sequential and mirror 
leptons (see Section 2.3.1). Using centre-of-mass energies up to 209 GeV, the L3 
experiment [322] was able to exclude new leptons with masses below around 103 
GeV to 95% CL. Varying the exotic scenario changes the limit by less than 1%. 

A number of searches have been made by the LEP experiments for stable SMPs 
within SUSY scenarios. ALEPH [325] has obtained lower-mass limits for stable 
coloured sparticles through an analysis of LEP-1 and LEP-2 data. In Ref. [325] 
MSSM scenarios with a stable gluino or squark and R-parity conservation was 
assumed (see Sections 2.1.1- 2.1.4). For a stable gluino hypothesis, the process 
e+e" qqgg with LEP-1 data was considered. In this scenario a radiated gluon 
splits into two stable gluinos. The hypotheses that i?-hadrons would be formed as 
either as neutral or charged states were both considered. In the former case, R- 
hadron jets were identified via a discrepancy between their measured hadronic en- 
ergy loss and the measurement of the total momentum of charged particles within 
the jet. A range of topological event shape variables were also used to suppress 
background. In the latter case ^ measurements, together with topological vari- 
ables, were used to identify i?-hadron events. A lower-mass limit for gluinos of 
27 GeV was evaluated at 95% CL with this work. A similar analysis by DEL- 
PHI [326] gave a lower limit of 18 GeV. Differences in the limits arise due to 
several factors. ALEPH used around 4 million hadronic Z° events compared with 
1.6 million events used by DELPHI. Furthermore, different topological variables 
were considered and different analysis techniques used to suppress background. In 
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addition, ALEPH employed a model of i?-hadron scattering in material which is 
based on the geometric cross-section, whilst DELPHI used the first BCG ansatz 
(see Section 5.2). The gluino mass limits of ALEPH and DELPHI were extracted 
assuming a direct production mechanism and are applicable to any MSSM scenario 
containing long-lived gluinos. Thus, these works can be considered to belong to the 
second class of searches. 

The ALEPH work also sought evidence of stable stops and sbottoms via the reac- 
tion e+e" qq, using ^ measurements. Since the partners of the left and right 
handed states of the top (or bottom) mix to form mass eigenstates, the lower mass 
limits depend on the choice of mixing angle [12, 13]. In this study, a stop (sbottom) 
mixing angle of 56° (68°) was used. At these values the stop and sbottom cou- 
plings to the Z vanish. This allowed the extraction of conservative mass limits: 95 
GeV (stop) and 92 GeV (sbottom), at 95% CL. As for the stable gluino work, the 
squark mass limits are relevant for any MSSM scenario accommodating long-lived 
squarks. 

Similarly, stable heavy slepton limits which are valid for all MSSM scenarios in 
which the slepton is stable have also been evaluated. Direct pair-production mech- 
anisms e+e" and e+e~ were assumed in a DELPHI search [294]. 
Within the MSSM, the cross section for each of these processes depends only on 
the relevant slepton mass. Lower-mass limits for the superpartners of the left (right) 
muon and tau of 98 (97) GeV were obtained at 95% CL. Mass limits of a sim- 
ilar value were also obtained by ALEPH [327] and OPAL [328]. Mass limits on 
selectrons are very model dependent owing to an additional t-channel neutralino 
exchange process. 

Searches for sleptons are usually interpreted within a slepton NLSP or a sleptons 
co-NLSP GMSB scenario, as described in Section 2.1.2. Such class 3 results ex- 
clude regions in parameter space defined by quantities such as the slepton masses, 
the effective SUSY-breaking scale, and the gravitino mass [294,327,328]. Since 
stable slepton scenarios represent only a small part of the GMSB parameter space, 
we do not provide a detailed description of these results here. 

Another type of class 3 SUSY search looks for NLSP charginos which are nearly 
mass degenerate with the LSP neutralino [321,329-331] , as can occur in AMSB 
and gravity mediated SUSY-breaking scenarios (see Section 2.1.4). For mass dif- 
ferences of less than approximately the pion mass, the charginos are stable during 
their passage through the detector. A search by ALEPH [329] provided the most 
stringent lower-mass limit of 101 GeV for stable charginos at 95% CL. However, 
as for the GMSB searches, this result is dependent on the SUSY model parameters. 
Refs. [321, 329-331] provide a detailed description of the model parameter spaces 
excluded by collider studies. 
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7.1.2 Searches at lepton-hadron experiments 

There have been few searches for electrically charged SMPs in lepton-hadron scat- 
tering experiments. Those which have taken place are generic (class 1) searches. A 
search for free quarks was performed by the EMC experiment [332] in low-energy 
deep-inelastic /x-Beryllium scattering (DIS). This experiment was sensitive to val- 
ues of the photon virtuality, Q"^, and the invariant mass of the entire hadronic final 
state, W, which extended to around 100 GeV^ and 20 GeV, respectively. Sets of 
scintillator counters were used to gain sensitivity to the anomalously low ionisation 
energy loss expected from free quarks. The ratio of the upper limit on the single 
inclusive production cross section for SMPs to the total inelastic muon cross sec- 
tion was determined to be around 10^^. The experiment was sensitive to SMPs 
with charges ±ie and ±|e for masses up to 15 GeV and 9 GeV, respectively. The 
limits on the production cross-section were smaller than those obtained in neutrino- 
nucleon scattering [333-335]. 

There has been one recent dedicated search for SMPs produced in high-energy 
lepton-hadron scattering, which was made by the HI experiment [275]. This used 
a data sample corresponding to processes with an average W of 200 GeV, and 
in which the exchanged boson was quasi-real (photoproduction). In photoproduc- 
tion processes, the exchanged photon can be ascribed a partonic structure, and, in 
this picture, the HI study is therefore equivalent to the hadron-hadron searches de- 
scribed in 7.1.3. However, there are differences between final states produced by 
photoproduction and hadron-hadron interactions, for example, due to the contribu- 
tions of so-called direct and anomalous photoproduction processes [336]. Thus, the 
two environments are not identical, and this search is therefore complementary to 
searches performed at a hadron-hadron facility. 

The HI experiment used measurements of ^ in its drift chamber to search for 
SMPs depositing anomalously high amounts of ionisation energy. A minimum bias 
sample corresponding to 6 pb~^ of luminosity was used. At 95% CL an upper cross- 
section limit of 0.19 nb was extracted for the production of SMPs. No production 
model was assumed for this study and the mass range explored is determined by 
the capability of the HI tracking chamber. Since the apparatus was fully sensitive 
to SMPs with charge ±e in the region of transverse momentum and mass 0.2 < 
Pt/M < 0.7, this implies a sensitivity in mass of up to around 100 GeV, assuming 
an SMP pair-production mechanism. 



7.1.3 Searches at hadron-hadron experiments 

Early class 1 searches for fractionally charged particles produced in pp interactions 
at the ISR (y^ = 52 GeV) were able to set upper cross-section limits for particles 
possessing fractional charges ±|e,±|e, and ±|e, and masses less than around 20 
GeV [337] . A host of lower energy hadron-hadron experiments also performed such 



68 



searches, usually for masses less than around 10 GeV [1]. 

Searches for SMPs in pp collisions at the Tevatron have been made by the CDF 
collaboration [338-340]. These were based on relatively low-luminosity samples 
of 26 nb-i [338], 3.5 pb"^ [339] and 90 pb^^ [340]. By the end of the Tevatron 
Run-II data-taking period, the DO and CDF experiments should each have collected 
several fb^^ of luminosity and we anticipate that the above mentioned studies be 
updated. Here we consider the works in Refs. [339, 340] which supersede the early, 
very low-luminosity CDF study [338]. 

In Ref. [339] SMPs leaving a slow (0.25 < (3 < 0.65) muon-like track were 
sought. The analysis also relied on measurements of ionisation energy loss, muon- 
like tracks and time of flight information from the calorimeter. Model-independent 
upper limits on the production cross section (calculated to 95% CL) of between 
120 and 5 pb were derived for the pair production of fractionally charged fermions 
of masses between 50 and 500 GeV. The electric charges considered were ±|e,±e, 
and ±|e. 

Using theoretical models, these data were also used to extract mass limits under 
the assumption that the SMP possessed a specific colour charge (singlet, octet or 
decuplet). These, and the remaining search results described here belong to the 
second classification of searches. The mass limits for an SMP with a given colour 
and electric charge vary between 140 and 255 GeV at 95% CL. 

The third CDF work [340] did not employ time of flight information, and relied on 
searches for particles losing anomalous amounts of ionisation energy. This study 
considers the possibility of fourth-generation quarks (see Section 2.3.1), and was 
able to set lower-mass limits at 95% CL for exotic quarks with charges ±|e and 
±|e at 220 GeV and 190 GeV, respectively. Furthermore, this work considered the 
production of fourth generation quarks without hadronisation and charge exchange 
effects. The limits obtained (0.3-2 pb over a mass range 100-270 GeV) broadly 
correspond to a generic search for the pair production of particles with charge ±e, 
assuming a strong production mechanism. 

Stable gluinos have not been explicitly considered in the Tevatron works. How- 
ever, it has been argued [341] that existing Tevatron limits on stable charged par- 
ticles [340] and anomalous mono-jet production [342] imply lower limits on the 
gluino mass to be roughly 170-310 GeV; the lower (upper) value corresponds to 
an extreme scenario in which gluino i?-hadrons are produced solely as neutral 
(charged) states. 

A heavy slepton hypothesis which assumed a Drell-Yan like production mechanism 

was also studied in [340] and production cross-section limits were extracted. How- 
ever, the expected cross section in exotic scenarios such as GMSB is more than 
an order of magnitude below this level of sensitivity. Contrary to the situation for 
exotic quark searches, heavy lepton mass limits from SMP searches have, so far. 
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only been extracted by the LEP experiments. 



7.1.4 Summary and discussion of direct searches 

A selection of class 1 searches is summarised in Tab. 4. Since these studies focus 
on generic SMPs, for which the production mechanisms are a priori unknown, it 
is important that searches were conducted in each of the different collision envi- 
ronments. Evidence for particles with fractional charges ±|e, ±|e, ±e, ±|e, and 
±|e was sought. Plastic track detectors deployed at LEP extend the SMP charge 
sensitivity to around |240e|. The mass sensitivity at the various accelerators for 
the pair production of generic SMPs with charge ±e extended to 100 (HERA), 101 
(LEP), and 500 GeV (Tevatron). None of the above generic searches were sensitive 
to particles with charge much below ~ ||e|. Dedicated searches for millicharged 
particles with charges as low as ~ 10~^e have been made at accelerators and else- 
where albeit typically for particles with masses less than ~ 1 GeV [3,98]. These 
searches are thus beyond the scope of this report. It is relevant for the massive par- 
ticles considered here to note that there exist no direct experimental searches which 
are sensitive to particles with charge ^||e| and mass 1<M<1000 GeV, and that 
indirect astrophysical constraints tend only to exclude such SMPs with charges less 
than around llQ-^e] [99]. 

For the class 2 searches, lower mass limits were obtained when considering a spe- 
cific type of SMP and an assumed production cross section. A compilation of lower 
limits for masses of specific types of particles is shown in Fig. 21. The stable quark 
limits are taken from Ref. [340] which assumed a fourth-generation quark model. 
The stable squark and gluino mass limits from LEP [325] are ~ 90 and 27 GeV, re- 
spectively. Stable heavy squark or gluino hypotheses have not yet been considered 
by Tevatron experiments although it could naively be expected that, in a minimal 
SUSY scenario with squark or gluino LSPs, mass limits with Tevatron data would 
be similar to those already obtained for stable quarks. Stable stau and smuon mass 
limits of around 98 GeV have been obtained at LEP [294, 327, 328]. Non-SUSY 
heavy lepton mass limits have also been extracted at LEP [292,322,343]. The mass 
limit is largely insensitive to the type of new lepton, and we quote here the L3 mass 
limit of 103 GeV [322] for mirror leptons. A summary of selected searches made 
within classes 2 and 3 is given in Tab. 5. 



7.1.5 Indirect collider searches 

Indirect mass limits of up to around ^ for fourth generation quarks and leptons, 

squarks (assuming a non-vanishing coupling to the Z°), sleptons, charginos, and 
exotic particles with electric charges as low as ~ 0.2e can be inferred from mea- 
surements of the invisible width [1,99]. Each limit is generally applicable for 
a scenario in which SMPs are pair produced through the decay of a and not 
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Collision Experiment Mass range Charge (e) Cross-section CL (%) Ref. 



(GeV) (GeV) Umit (pb) 



1800 


PP 


CDF 


100-270 


±1 


0.3-2 


95 


[340] 


1800 


PP 


CDF 


50-500 


^3 


10-100 


95 


[339] 










±1 


5-50 
















5-70 






300 


ep 


HI 


< 100 


±1 


190 


95 


[275] 


130-209 


e+e~ 


OPAL 


45-102 


=^3 


0.005-0.02 


95 


[280,319] 








45-102 


±1 


0.005-0.03 












45-100 


±1 


0.005-0.02 












45-98 


±f 


0.005-0.02 






189 


e+e~ 


DELPHI 


68-93 


±1 


0.02-0.04 


95 


[292] 


130-183 


e+e~ 


DELPHI 


2-91 


±1 


0.05-0.3 


95 


[293] 


130-183 


e+e~ 


DELPHI 


2-91 


=^3 


0.04-0.6 


95 


[293] 


130-172 


e+e^ 


ALEPH 


45-86 


±1 


0.2-0.5 


95 


[276] 


91.2 


e+e" 


ALEPH 


5-45 


±1 


3 - 10 


90 


[279] 










±1 


1-12 












10-72 


±1* 


1.6-140 














=^3 


1.4-4 






91.2 


e+e~ 


OPAL 


3-45 


=^3 


0.2-1.0 


90 


[318] 










±2* 
=^3 


0.15-0.9 


95 












±1* 


0.15-3.0 














^3 


0.18-0.21 














±2* 


0.27-0.3 






91.2 


e+e~ 


Kinoshita et al. 


1-45 


<|240e| 


70 


95 


[286] 



Table 4 



A summary of selected direct searches for electrically charged SMPs belonging to class 1 . 
The searches are categorised according to centre-of-mass energy, colliding particles, and 
experiment. A range in cross-section limit is provided for an SMP with a given charge. 
The mass range corresponds to the region for which the upper cross-section limit is quoted. 
The limits were derived under the assumption of a pair-production mechanism, with the 
exception of those results marked with the symbol *. For these limits, a single SMP inclu- 
sive production mechanism was assumed. The confidence level to which the limits were 
extracted is also shown. 

recorded in a detector. Whilst the majority of the aforementioned particles would 
certainly be detected in the direct searches, it is, however, conceivable that frac- 
tionally charged particles may have been missed. LEP-1 limits for exclusive pair 
production processes were not made for SMPs with charges below ±|e. Further- 
more, in a free quark picture, these particles may possess a large interaction cross- 
section, and rapidly become stopped in detector material. In these scenarios, the 
indirect limits provide a valuable complement to the direct searches. 

Similarly, indirect mass limits can be used to support the stable gluino limits from 
direct searches, which are dependent on the phenomenological models used to de- 
scribe i?-hadron energy loss. These models possess theoretical uncertainties which 
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V3(GeV) 


Collisions 


Experiment 


Particle 


Mass limit (GeV) 


Ref. 


1.8 TeV 


PP 


CDF 


4*'' generation quark 


rrit' >220, My >190 


[340] 


130-209 


e+e" 


OPAL 


Heavy leptons 


m;/ > 102 


[280] 


133-208 


e+e" 


L3 


Heavy leptons 


mil > 103 


[324] 


91.2-209 


e+e~ 


ALEPH 


Squark 


rrif > 95, mj > 92 
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L3 
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Table 5 



A summary of selected classes 2 and 3 searches which have been made for electrically 
charged SMPs. The searches are described according to the collision centre-of-mass energy, 
colliding particles, experiment, the type of SMP, and the resultant mass limit. The squark, 
gluino, and slepton searches are valid for any MSSM scenario in which these particles 
are long-lived. The squark and slepton limits correspond to values of mixing angles which 
predict the lowest possible cross-sections. The chargino limits are model-dependent and 
were extracted in AMSB and AMSB-like scenarios. All limits were extracted at 95% CL. 

are not trivial to estimate. Indirect lower limits on the masses of colour octet par- 
ticles of around 6 GeV can be set from the partial width of any purely hadronic 
exotic contribution (negative or positive) to Z-decays [344]. 

Mass constraints on gluinos have also been set by studying the QCD colour factors 
from 4-jet angular correlations and the differential 2-jet rate in Z-decays [345]. In 
addition, the modification in the running of ag as a consequence of the existence 
of stable coloured particles has been used to set limits on the mass of these parti- 
cles [346]. These lower mass limits also extend to around 6 GeV. 



7.2 Searches for magnetically charged particles 

A body of searches has been performed for magnetic monopoles produced at accel- 
erators, in cosmic rays, and bound in matter [1]. Although no reproducible evidence 
for the existence of monopoles has been found, several apparent observations have 
intermittently caused great excitement in the field. In 1973 a particle consistent 
with being a monopole with charge around 2g£) was recorded at a cosmic ray ex- 
periment [347]. While the event remains to be fully understood, the authors later 
established that a monopole could not be responsible for the observation [348]. 
In the 1980's isolated observations [349, 350] of cosmic ray monopole candidates 
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Fig. 21. Lower mass limits for stable gluinos [325], squarks [325], fourth-generation 
quarks [340], and heavy leptons [322]. The lengths of the vertical Unes have no signifi- 
cance. 

were made using SQUID detectors. However, instrumental effects could not be 
ruled as being responsible for the signals and the observations are inconsistent with 
monopole flux limits set at other cosmic ray facilities [351]. 

During the 1990's monopole searches remained a routine part of the experimental 
programs of collider experiments. This section discusses searches for monopoles 
directly produced at electron-positron annihilation, hadron-hadron and lepton-hadron 
facilities. Direct searches include, for example, searches for monopoles leaving 
tracks in plastic track detectors, and searches for monopoles trapped in detector 
material. As will be shown, most searches are concerned with the possibility of 
monopoles possessing the Dirac charge. However, even in the absence of theoret- 
ical motivation, it is prudent to also consider scenarios with particles possessing 
as wide a charge-range as is experimentally possible. We therefore describe the re- 
gions in magnetic charge and mass which have been explored. The discussion will 
be followed by a brief discussion of indirect monopole searches. 



7.2.1 Searches at e experiments 

A variety of methods have been used to search for monopoles at e+e^ experi- 
ments. The most commonly used approach, the deployment of plastic track de- 
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lectors around the beam interaction point, has been used in searches at the PE- 
TRA [352], PEP [353,354], KEK [355,356] and LEP-1 [253,286] facilities. The 
most competitive exclusion limit for Dirac monopoles was extracted with passive 
detectors deployed at the OPAL intersection point at LEP-1 [253]. No monopoles 
were observed and an upper cross-section limit for the production of Dirac monopoles 
of 3 nb was computed at 95% CL. Since the acceptance depends on the monopole 
mass and charge, the maximum value of acceptance was used in this extraction. It 
was also assumed that monopoles would be exclusively produced in pairs and dis- 
tributed isotropically in phase space. The experiment was sensitive to monopoles 
with magnetic charge between 0.9gn and S.6gD and for masses up to 45 GeV. 

An earlier experiment at LEP by the MODAL [357] group used CR-39 plastic sheet 
detectors at the 15 LEP interaction point [286], as mentioned in Section 7.1.1. As in 
the case of [253] this experiment was sensitive to monopole masses up to around 45 
GeV albeit with poorer sensitivity to the monopole production cross section (a < 
0.1 fib at 95% CL for Dirac monopoles). However, the experiment was sensitive to 
lower values of magnetic charge (0.1 Lower-energy plastic track experiments at 
e+e" facilities have also been made albeit with a lower mass sensitivity [352-356]. 

Searches for monopoles leaving parabolic tracks have been made [298, 358] (see 
Section 6.4.1). These assumed an exclusive pair-production mechanism. Together, 
these works gave sensitivity to values of the magnetic charge between 0.03 — O.lgf^ 
for monopole masses below 17 GeV. 

7.2.2 Searches in lepton-hadron collisions 

There has been only one search for monopole production in lepton-hadron scat- 
tering. A search was made for monopoles which could have been stopped in the 
aluminium beampipe of the HI experiment at HERA [254]. The beampipe was 
exposed to an integrated luminosity of 62 pb"^ The SQUID apparatus, which was 
used in this experiment is shown in Fig. 19. To calculate the acceptance in this study 
it was assumed that monopoles were produced in pairs via photon-photon fusion 
processes. Limits were extracted for two different assumed reactions: ep eMMp 
and ep — > eMMX. In the former (latter) case the monopoles were treated as 
spin (i) objects. No evidence was found for the presence of trapped magnetic 
monopoles. This work set mass-dependent upper limits on the production cross 
section of monopoles with charges gD,'2gD,^gD, and dgoior masses up to around 
140 GeV. The cross-section limits (calculated at 95% CL) vary between around 0.2 
and 100 pb depending on the charge, mass and assumed monopole spin. 

7.2.3 Searches at hadron-hadron experiments 

A recent study at CDF [252] obtained upper cross-section limits for monopoles pos- 
sessing the Dirac charge for masses up to 900 GeV. This work used 35.7 pb~^ of 
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integrated luminosity and employed a specially built highly ionising particle trigger 
requiring large light pulses at both ends of a TOF scintillator bar. In addition, track 
reconstruction software was optimised to search for characteristic parabolic trajec- 
tories. A Drell-Yan-like pair-production mechanism, as described in Section 4.6, 
was used to simulate the kinematic properties of monopoles. Using this model, 
cross sections of greater than 0.2 pb (at 95% CL) were ruled out for masses be- 
tween 200 and 700 GeV. Using the same Drell-Yan formalism as a cross-section 
prediction, this can be interpreted as a lower mass limit of 360 GeV. 

Another recent search at the Tevatron involved the use of the SQUID technique by 
the E882 experiment to detect monopoles bound in detector material at the DO and 
CDF experiments [359,360]. In these works three sets of samples were taken from 
the DO and CDF detectors during their upgrades for Tevatron Run-2 data-taking: 
the Be beampipe and Al extension cylinders from DO, Pb from the CDF Forward 
Electromagnetic Calorimeters, and one half of an Al cylinder supporting the CDF 
Central Tracking Chamber. These detector pieces were exposed to integrated lumi- 
nosities of 172 ± 8 pb-i (DO) and 180 ± 9 pb'^ (CDF). Three different assumptions 
were made concerning the angular distribution of monopoles in the monopole- 
antimonopole centre-of-mass system: oc constant, oc 1 + cos'^6, and 
cx 1 — cos'^9. Upper limits on the production cross section of monopoles with 
charge gd, 2gd, Sqd and Gqd were found to be 0.6,0.2, 0.07, and 0.02 pb, respec- 
tively, at 90% CL, when taking the isotropic case. Using the modified Drell-Yan 
formalism as a cross-section prediction, these correspond to mass limits of 265, 
355, 410, and 375 GeV, respectively. Varying the angular distributions according 
to the three assumed forms described above leads to a spread of around 10 GeV in 
these mass values. 



Taken together, the works [359, 360] and [252] represent the most comprehensive 
limits for the production of monopoles in hadron-hadron collisions. The searches 
each employ complementary methods, and together are sensitive to the production 
of monopoles with charges between go and Gqd and masses up to around 900 GeV. 
Nevertheless, since both CDF and DO are expected to accumulate several fb~^ of 
luminosity at the completion of the Tevatron program, there remains sensitivity to 
search for monopoles with a cross section around an order-of-magnitude lower than 
that which has already been excluded. In addition, a further search could address 
the possibility of the production of monopoles with a charge less than gr,. 

Other searches [282-284] at the Tevatron based on substantially smaller samples 
of luminosity have used plastic track detectors placed around the interaction point. 
Again, no evidence for monopole production was found and the exclusion limits 
on cross sections are typically several thousand times larger, and the mass limits 
around three times smaller than those obtained in [252, 359, 360]. 

Lower energy hadron-hadron experiments have employed a variety of search tech- 
niques including plastic track detectors [361,362] and searches for trapped monopoles 
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[306-308]. The latter searches used external magnetic fields to accelerate stopped 
monopoles and were sensitive to a very wide charge range (0.03-24g'£,) albeit for 
masses less than 30 GeV. 



7.2.4 Summary and discussion of direct searches 

In Tab. 6, the upper cross-section limits from monopole searches and their charge 
and mass sensitivities are given. It is difficult to compare the mass and cross-section 
exclusion limits extracted in the various searches because of the theoretical uncer- 
tainties on monopole production. What can be reliably calculated is the region of 
finite acceptance for a monopole with a specific charge and mass, and then, with 
some model dependence, the derived upper limit on the production cross section. 
Mass limits, extracted under the assumption of a known production cross section, 
are the least reliable information reported by experiments. 

As motivated by Dirac's argument (see Section 2.4.1), a large number of experi- 
ments have searched for monopoles possessing the Dirac charge. Fig. 22 shows the 
upper cross-section limits, which have been obtained. The limits are shown versus 
one half of the centre-of-mass energy of the experiments, which is indicative of 
the largest mass value to which an experiment is sensitive, currently 900 GeV at 
the Tevatron [360] . There also exist theoretical arguments for considering parti- 
cles with charges up to Gqd- To reach these values it is usually necessary to look 
for stopped monopoles owing to the vastly increased energy loss associated with 
high-charge monopoles. Searches for monopoles with magnetic charges up to Ggo 
have so far only taken place at two high energy colliders: HERA [254] and the 
Tevatron [360]. 

Sensitivity to low magnetic charges (^^^/d) was provided by the HI [254] and 
MODAL [286] searches. Low-energy e+e^ experiments searching for parabolic 
tracks [298, 358] and hadron-hadron searches in which external magnetic fields 
were used to accelerate trapped monopoles [306-308], provided sensitivity to charges 
as low as O.OSgD- However, these searches were only sensitive to monopoles of 
masses less than around 30 GeV. 

There have been no dedicated searches for monopoles with charges less than O.OSgn- 
Searches for electrically charged particles with anomalous ionisation energy loss 
(Section 7.1) could be sensitive to monopoles in this low charge region. However, 
this depends on the dynamic range and calibration of the read-out electronics of 
the tracking devices, which were used in these searches. In addition, measurements 
of charged particles with high transverse momentum could also provide sensitivity 
to low charge monopoles. Since a monopole with a large mass (> 100 GeV) and 
low charge would not leave an identifiable parabolic track, the signature would be a 
straight track consistent with a very high-momentum electrically charged particle. 
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Table 6 



A summary of selected direct searches for monopoles. The searches are categorised accord- 
ing to the collision centre-of-mass energy, the colliding beams, the experiment, the quoted 
acceptance regions in charge and mass, and the cross-section Umits for specific values of 
magnetic charges. The symbol "-" denotes that an experiment did not explicitly quote an 
acceptance region. The cross section limit from Ref. [252] is quoted for monopole masses 
between around 200 and 700 GeV, a region in which the limit is roughly constant. The 
cross section limits from Ref. [254] are quoted for a monopole mass of 100 GeV, which is 
approximately one half of the average photon-proton centre-of-mass energy at the experi- 
ment. The cross section limits for works [298] and [358] are quoted for mass values of 8 
and 2.5 GeV, respectively. At these values cross-section limits are available for the different 
charge hypotheses considered. 

7.2.5 Indirect searches 

IVIass sensitivity at colliders can be increased via indirect searches, which rely on 
processes mediated by virtual monopole loops, as shown in Fig. 9 in Section 4.6. 
Such processes can lead to multi-photon final-states arising from electromagnetic 
radiation from monopoles. The DO experiment were able to exclude monopole 
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Fig. 22. Summary of upper limits on the monopole production cross section at a range of 
experiments at different centre-of-mass energies [252-254, 308, 352, 354, 356, 359-361]. 

masses between 610 and 1580 GeV [258], while L3 reports a lower mass limit 
of 510 GeV [257] using this approach. However, it should be stressed that these 
cross-section calculations employ perturbative techniques [255, 256] with uncer- 
tainties which are difficult to estimate and the derived cross-section limits have 
been criticised for this [259]. 

An indirect lower limit on the mass of Dirac monopoles of 120 GeV has also 
been inferred from measurements of the anomalous magnetic moment of the muon 
[363]. However, it should again be emphasised that this result relies on perturbative 
techniques and that the result is therefore questionable. 



8 SMPs at the LHC 



The LHC will allow searches for SMPs with masses up to several TeV, repre- 
senting an order of magnitude increase in mass sensitivity compared with earlier 
colliders. As explained in Section 2, well motivated extensions of the SM predict 
the production of a range of SMPs at the LHC. However, it is important that the 
LHC experiments perform searches for SMPs which are motivated not only by 
what is currently theoretically favoured but also by what is experimentally feasible. 
Therefore, in this section we discuss the sensitivity of the two high-luminosity and 
multi-purpose experiments, ATLAS [364,365] and CMS [366,367], to SMPs with a 
variety of properties. In addition, we also briefly discuss the planned MOEDAL ex- 
periment [368,369], which aims to deploy plastic track detectors around the LHCb 
vertex detector. 
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We use the findings of a number of studies which explore the feasibility of discover- 
ing SMPs at the LHC. These have been made both by interested phenomenologists 
and the experiments themselves (see for example [243,341,370-380]). The for- 
mer works tend to apply simple ansatzes regarding detector performance whereas 
the latter use sophisticated detector simulation packages. While the latter studies 
provide a more realistic description of the detector response it is pertinent to note 
that the detector simulation programs are still being optimised and have yet to be 
confronted and verified with collision data. Nevertheless, it is still possible to use 
the feasibility studies to draw a number of conclusions regarding the discovery po- 
tential for various SMPs and to highlight SMP-specific experimental challenges. 



8.1 Detector sensitivity 
8.1.1 Timing issues 

One of the principal constraints on the design of a collider experiment and its ability 
to measure slow-moving SMPs is the bunch crossing time. The LHC crossing time 
of 25 ns, will be considerably shorter than that at LEP (25 //s), the Tevatron (396 ns) 
or HERA (96 ns). Therefore, for an SMP at the LHC to be detected or triggered in 
a certain detector system and be associated to the correct bunch crossing, it should 
arrive at most 25 ns after the default arrival time of a particle travelling at the speed 
of light [377, 38 1]. Later arrival would imply triggering or detection within the next 
crossing time window. The large size of the ATLAS and CMS detectors (the central 
ATLAS and CMS muon chambers extend to 10 and 7 m, respectively) ensures that 
this will be an important source of inefficiency in detecting SMPs. For example, 
it is only possible to reconstruct the track of a slowly moving SMP in the ATLAS 
central muon chambers within the correct bunch crossing window if /3>0.5 [377]. 
Even if the SMP travels within the appropriate timing window, additional problems 
may arise from its slowness. The sampling time and reconstruction software of each 
sub-detector is optimised assuming particles travelling at luminal speed. Hence, the 
quality of the read-out signal or reconstructed track or cluster may be degraded for 
an SMP, especially for sub-systems far away from the interaction point. Detector 
simulations so far suggest that it will still be possible to trigger and measure slowly 
moving particles at ATLAS and CMS [372, 377]. However, this is an area which 
must continue to be studied as the simulation programs are further developed and 
the detectors better understood. 



8.1.2 Important sub-detectors 

The principal detector systems needed for SMP searches are the tracking systems. 
ATLAS and CMS both contain inner tracking systems capable of precision mea- 
surements of the momenta of tracks of charged SMPs [364,366]. These devices can 
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potentially be used to search for the 'classic' signature of SMPs losing anomalous 
amounts of ionisation energy in an inner detector, albeit probably with poorer pre- 
cision than has been obtained in previous experiments. The ATLAS inner tracking 
system comprises a semiconductor tracker system and a transition radiation tracker 
(TRT). The precision with which ^ can be exploited is limited since the TRT will 
only record hits satisfying high and low ionisation thresholds. The total number 
of high/low threshold hits and the time in which a channel is read out over a spe- 
cific threshold can be used in SMP searches [377]. The CMS inner tracking system 
consists of silicon pixel and silicon strip detectors. As described in Section 5, the 
resolution in ^ of such devices is poorer than with gaseous detectors. Work is 
ongoing to estimate the sensitivity of the CMS tracking system to SMPs [277]. 

Arguably the most important detector component in an SMP search at the LHC is 
the muon system since most proposed SMPs leave a signature of a slow-moving 
penetrating particle (see Sec. 6). The ATLAS [364] and CMS [366] experiments 
employ muon chambers which will be able to trigger on SMPs efficiently. They 
also offer considerably greater precision in track reconstruction than has been pre- 
viously available at current colliders experiments [382-385]. The excellent mo- 
mentum resolution implies that charge misidentification will be typically less than 
1% [364], and that i?-hadrons undergoing the charge exchange reactions described 
in Section 5.2 can therefore be identified [378]. In addition to providing the recon- 
struction of the track of an SMP, the muon systems will be able to provide time of 
flight measurements with an accuracy of around Ins. The time of flight is likely to 
be one of the most powerful discriminants between SMPs and muons [372, 377]. 
Furthermore, in the event of a discovery, it will allow the reconstruction of the 
mass of an SMP. It has been shown that a mass resolution of around 20 GeV is 
possible for an i?-hadron of mass 500 GeV at ATLAS [379] using this method. An 
understanding of the SMP mass resolution is important in both probing the mass 
hierarchy of particles in any exotic scenario and also in establishing the existence 
of different types of SMPs with similar masses. In addition to the time of flight 
information, the muon systems employ gaseous detectors which can provide ^ 
measurements. These may also be used as discriminants in an SMP search and in 
measurements of SMP properties. 

Since direct searches for most species of SMPs are based on charged particle sig- 
natures, calorimeters are less crucial, though still important, detector components. 
The ATLAS and CMS calorimeters can be employed in a number of ways in SMP 
searches. Measurements of electromagnetic and hadronic energy loss can distin- 
guish between hadronic and leptonic SMP species. Furthermore, the electromag- 
netic calorimeters can also provide supplementary ^ information. Calorimeter 
triggers based on jet multiplicity or global event properties, such as visible energy, 
can also be used to select events containing SMPs [377]. These same observables 
could be used in offline analyses for SMPs predicted to occur in specific event 
topologies. 
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A complementary way to detect SMPs at the LHC is via passive detectors. The 
MOEDAL collaboration [368,369] proposes to deploy a plastic track-etch detector 
comprising three layers of CR39/lexan plastic around the LHCb vertex detector. 
The major systematic uncertainty in any search with MOEDAL will be background 
tracks from the radiation environment of the LHC, and this is presently under study. 

8.2 Specific types of SMPs at the LHC 

8.2.1 Expected rates 

Using Pythia we have calculated the expected yields for different types of SMPs 
at the LHC, following the accumulation of 10fb~^. This corresponds to one year 
of low-luminosity operation at the LHC. Table 7 shows the expected yield of pair- 
production processes giving rise to different types of SMPs: gluinos, stops, fourth- 
generation quarks with charge ±|e, heavy leptons and superpartners of left and 
right-handed leptons. Since predictions of this type are necessarily scenario-dependent, 
we have made conservative calculations assuming that the exotic particle is pro- 
duced via SM couplings and that the influence of other exotic particles is minimal. 
The predictions do not include contributions from the cascade decays of other ex- 
otic particles. 

Studies indicate that the acceptance of the ATLAS and CMS detectors following 
SMP selection cuts typically vary between 5% and 80% depending on the mass 
and type of SMP [370-375, 377, 378]. A detailed simulation for each type of SMP 
of detector inefficiencies and systematic uncertainties, such as the effects of nuclear 
scattering and delayed reconstruction, is ongoing and therefore beyond the scope 
of this paper. Furthermore, many systematic uncertainties, particularly for searches 
for hadronic SMPs, remain difficult to quantify even with detector simulations. This 
is discussed below, where we describe experimental issues related to specific types 
of SMPs which could be sought at the LHC. 
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Table 7 

Pythia predictions of the total number of SMP pairs expected at the LHC following the 
accumulation of 10fb~^ of integrated luminosity. 
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8.2.2 Generic SMPs 



Generic SMP's which possess electric charge and penetrate through the whole de- 
tector will leave a clear signature of a high-p^, slow, muon-like track. Since no 
assumption, other than non-zero electric charge, is made regarding the properties 
of a generic SMP it is not possible to calculate its expected production cross sec- 
tion. However, background is expected to be highly suppressed when making use 
of the time of flight technique [377]. Thus, penetrating SMPs with appropriate pro- 
duction cross sections would be easily observed. A null observation would also 
allow stringent upper limits on these cross sections. For SMPs which stop in the 
calorimeter and do not traverse the muon system, searches could rely on observ- 
ing an excess of high-p^ tracks in the inner detectors. In conjunction with this, 
observables such as a charged particle leaving a signature of anomalous ionisation 
in the tracking chambers together with a characteristic energy deposition profile in 
the calorimeters could also be used. However, as described earlier in this section, 
the precision with which ionisation energy loss can be measured is likely to be 
poorer than obtained in earlier searches. This could make it difficult in particular to 
identify fractionally charged objects such as for example free quarks. 

The planned MOEDAL experiment would be able to complement the above search 
strategies by allowing searches for electrically charged SMPs with ^ > 5, which 
implies a maximum electric charge sensitivity of up to around 220e [386]. The 
present upper limit on the SMP electric charge to which ATLAS and CMS will be 
sensitive is still unclear. However, it is likely to be substantially below the reach of 
MOEDAL. 

While it can be expected that the ATLAS and CMS tracking systems will be able 
to reconstruct particles with charges down to at least ~ 0.5e, studies have not been 
performed to establish this and to estimate the reconstruction efficiency. Measure- 
ments of SMPs with extremely low charge (millicharged-charged particles) would 
not be directly recorded in the LHC tracking chambers. 



8.2.3 Heavy leptons and sleptons 

The possibility of stable heavy leptons and sleptons at the LHC hs been considered 
in Refs. [371-376,376,380]. As these studies illustrated, the detection of these 
particles presents the fewest detector difficulties of all of the SMPs considered, al- 
though their direct cross section would be suppressed in relation to coloured SMPs. 
Heavy leptons and sleptons would manifest themselves as delayed, muon-like par- 
ticles leaving little energy deposition in the hadronic calorimeters. Since these par- 
ticles would not necessarily be associated with a jet, isolation criteria could also be 
used to identify them. In an extremely minimal scenario in which only the proper- 
ties of the slepton or lepton would play a role in the production cross section, mod- 
elling these process is relatively straightforward. However, the presence of other 
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exotic particles complicates the picture. In SUSY scenarios accommodating stable 
sleptons, e.g. GMSB, stable sleptons would be expected to be dominantly produced 
via the cascade decays of copiously produced squarks and gluinos at the LHC. Stud- 
ies of stable sleptons have focused on their detection and the exploration of the pa- 
rameter space of GMSB and GMSB-like models [371-376, 380]. For conservative 
estimates of rates, Tab. 7 presents the expected yields of the stable superpartners 
of left and right handed leptons, which are produced directly and not via cascade 
decays. The latter quantities depend on masses and couplings of many SUSY spar- 
ticles so are extremely model-dependent. Making use of the time of flight tech- 
nique (see Section 6.3), the discovery of leptons or sleptons possessing masses up 
to several hundred GeV should be possible in early LHC running with the ATLAS 
detector [371,372,374,375]. In the optimistic situation of an SMP-discovery, the 
angular distribution of the pair-produced SMPs could allow discrimination between 
a lepton or slepton hypothesis, as outlined in Ref. [376]. 



8.2.4 Gluinos 

Gluinos are produced through strong interactions and would thus be copiously pro- 
duced at the LHC, as is shown in Tab. 7. A number of studies exploring the pos- 
sibility of gluinos at the LHC have taken place [243, 341, 377, 378], which were 
mainly inspired by the recent split SUSY model. These studies indicate that stable 
gluinos possessing masses up to 1 TeV could be discovered in early LHC running. 
However, there are special issues associated with hadronic SMPs which must be 
addressed in a stable gluino search. For example, poorly understood nuclear inter- 
actions could potentially inhibit a discovery. It is possible, though extremely un- 
likely, that charged gluino i?-hadrons could dominantly be stopped or convert into 
neutral states in the calorimeter and thus escape detection in the muon chambers. 
Furthermore, even if the i?-hadrons remain charged then events could be produced 
containing tracks with different signs of charge in the inner and muon systems (so- 
called 'flippers'), possibly challenging the track reconstruction software. Current 
simulations suggest that the ATLAS reconstruction programs can deal with this sit- 
uation, and it may in fact be used as a means of searching for /?-hadrons [378]. 
However, in the absence of collision data to confirm this, a conservative approach 
would be to consider the most straightforward signature in which the i?-hadron 
maintains the same non-zero charge in the inner and muon chambers. The string 
fragmentation model with default parameters (Section 4.3) and the geometric scat- 
tering model (Section 5.2.2) predict that the proportion of gluino i?-hadrons pro- 
duced at the LHC which possess the same non-zero electric charge at the production 
vertex and after their passage through a typical calorimeter is around 10%. In this 
situation, gluino i?-hadrons could be found with masses up to at least 1 TeV in 
early LHC running when using measurements of muon-like tracks and topologi- 
cal variables, and significantly beyond when also making use of the time of flight 
technique [377]. However, it should be noted that calculations of nuclear scattering 
processes for i?-hadrons are extremely uncertain and any search must take this into 
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account. 



Another pessimistic scenario can be envisaged if the fraction of gluino-gluon states 
formed in the hadronisation process is around unity, leading to no track in the inner 
detector. In this unlikely, but not inconceivable scenario, neutral i?-hadrons which 
convert into charged states via nuclear interactions may be visible as tracks in the 
muon systems. However, as for the charge 'flippers' this signature may again chal- 
lenge reconstruction software and the event filtered away with background removal 
algorithms. 

The most experimentally difficult scenario occurs in the case of gluino _R-hadrons 
being produced as neutral particles and remaining neutral during their passage 
through the detector. Here, a signature of i?-hadron production could be manifest 
through an anomalously high jet rate [341]. 

8.2.5 Stable s quarks and quarks 

Stable squarks and quarks at the LHC have been postulated, as described in Sec- 
tion 2 although no detector simulations of these have been published. Searches for 
stable heavy colour-triplet states would benefit from the strong production cross 
section and exploit the same techniques as for the stable gluino searches described 
above. Similarly, they would also be sensitive to uncertain nuclear interactions 
which could potentially reduce their discovery potential. The nuclear interactions 
of stable colour triplets have received far less attention than gluinos and searches 
would benefit from more phenomenology in this area. One model which treats both 
colour-octet and triplet states [235] predicts comparable energy loss and charge ex- 
change rates for both types of particles. Should this model reproduce the data, it can 
therefore be expected that a similar sensitivity is obtainable as for stable gluinos, 
i.e. stable squark and quarks with masses up to at least 1 TeV could be discovered at 
early stages of the running of the LHC. Again, the angular distribution of the pair- 
produced SMPs could allow discrimination between a quark or squark hypothesis 
in case of a discovery. 

8.2.6 Monopoles 

Magnetic monopoles at the LHC will leave a number of striking signatures. How- 
ever, little work has been done estimating the sensitivity of ATLAS and CMS to 
these particles. Those studies which do exist concern indirect searches in which 
magnetic monopoles are produced as internal loops [365]. Furthermore, as dis- 
cussed in Section 4.6, calculations of these processes suffer from unquantifiable 
uncertainties [255,256]. Even in the absence of detailed studies of direct magnetic 
monopole production it is possible to make several remarks concerning their de- 
tection. As described in Section 5.1.2, a Dirac monopole will typically lose several 
thousand times as much ionisation energy as a MIR A detailed ^ calibration is 
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therefore not necessary to observe them. However, electronic saturation effects due 
to the enormous ionisation energy must be carefully studied in any search. Fur- 
thermore, customised track-finding algorithms may have to be written to account 
for the parabolic track trajectory followed by a magnetic monopole in a magnetic 
field. It is also important to employ GEANT simulations to calculate the stopping 
of magnetic monopoles in detector material and the energy deposition of profiles of 
magnetic monopoles which progress to the calorimeters. A comprehensive search 
must consider the full mass and magnetic charge range which is experimentally 
available. 

The MOEDAL experiment bypasses the problems described above of detecting 
magnetic monopoles with active detectors by using the plastic track technique. Cur- 
rent acceptance calculations indicate that MOEDAL will be sensitive to monopoles 
with masses up to around 6 TeV and charges up to around 3gD [386]. MOEDAL 
may offer the most promising method of hunting magnetic monopoles at the LHC. 
However, this depends on the magnetic monopole production cross section since 
MOEDAL would be exposed to integrated luminosity around 100 times lower than 
that collected by ATLAS and CMS. 



9 Summary 

Searches for Stable Massive Particles have been made at colliders for several decades. 
So far, no particles beyond those accommodated within the Standard Model have 
been observed. The searches are motivated by a number of theories which address 
key issues in modem physics. A review has been given of theoretical scenarios 
predicting SMPs, the phenomenology needed to model their production at collid- 
ers, the experimental techniques used to find SMPs and the searches which have 
been made to date. The interplay between collider searches and open cosmological 
questions has also been addressed. 

We look forward to the extension of these searches in the coming years at the LHC 
and at cosmic ray facilities. The discovery of an SMP would change our view of 
particle physics. 
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